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Equations differentielles 



1 Introduction 



Les oscillations lineaires des systemes a un degre de liberte sont regies par des equations differentielles du second 
ordre a coefficients constants. Une equation differentielle lineaire du second ordre a coefficients constants est 
une relation entre une variable dependante y et ses derivees premiere et seconde par rapport a une variable 
independante t, qui peut s’ecrire sous la forme : 



d 2 y 
dt 2 



+ + w o v — 



Les coefficients S et u>o sont des constantes reelles positives. Dans les problemes de vibration le parametre t 
represente le temps et on note par convention: 



dy 

dt 

d 2 y 

dt 2 

D’ou l’ecriture de l’equation differentielle : 



y 

V 



y + 25y + uly = A(t) 



2 Equation homogene 

L’equation differentielle est elite homogene si A(t) = 0: 

y + 25y + u 2 0 y = 0 

On recherche des solutions qui sont de la forme y(t) = e st . Dans ce cas l’equation differentielle s’ecrit : 

[s 2 + 26 s + Wq] e st = 0 
Cette equation doit etre vraie quel que soit t, ce qui implique : 

s 2 + 25 s + ujq = 0 

On obtient ainsi une equation du second degre en s dite equation caracteristique. Les racines de cette 
equation caracteristique sont : 

si = -S + \Juq- 5 2 




La solution de l’equation differentielle s’ecrit alors : 

y{t) = A ie Slt + A 2 e S2t 

A-[ et ff -2 sont deux constantes d’integration cpie Ton peut determiner a partir des conditions initiales: 

y (t = 0) = y 0 
y(t = 0) = y 0 
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Sachant que S et wo sont des nombres reels positifs, si et S2 sont negatives ou complexes avec une partie 
reelle negative. La nature des solutions Si et S2 de l’equation caracteristique depend de la valeur relative de <5 
par rapport a Wo ■ Ainsi trois cas sont a envisager : 

2.1 Regime fortement amorti ( <5 > oj 0 ) 

Dans ce cas, les deux racines Si et S2 sont reelles et negatives. L’ecriture des conditions initiales donne un 
systeme de deux equations : 



Ai + A2 — yo 

S1A1 + s 2 A 2 = y 0 

dont les solutions sont les constantes d’integration Ai et A2: 



S 2 yo - Vo 
S2 ~ si 
yo - Siyo 
S2 - Si 

y[t) est une fonction decroissant sans oscillations vers zero lorsque t augmente. La forme exacte de y(t) 
depend des valeurs de Ai et A2 qui sont determinees par les conditions initiales. Pour les conditions initiales 
particulieres suivantes :yo ^ 0; yo = 0, la figure suivante represente les variations de y au cours du temps 



Ai = 
A2 




2.2 Regime critique ( S = ojo ) 

L’equation caracteristique possede une racine double : 



Si = S2 = -6 

La solution generale de l’equation differentielle est de la forme : 

y(t) = (Ai + A 2 t) e~ St 

Les constantes d’integration Ai et A2 sont determinees a partir des conditions initiales et valent : 
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M = vo 
Ao = yo + <5 yo 

y(t) et une fonction decroissant vers zero quand t augmente. II est aise de verifier que cette situation 
correspond a la decroissance la plus rapide de y(t). Ce cas correspond au regime critique. Pour le cas particulier 
de conditions initiales :yo ^ 0; yo = 0, la solution est : 

V (t) = yo e~ st (1 + 5t) 

La figure ci-dessous represente les variations de y(t). 




2.3 Regime pseudo-periodique ( 8 < ojo ) 

Dans ce cas Si et s 2 sont complexes conjugues : 



si = -5 + j\Jul~ 5 2 
s 2 = -5- j\Jwl -5 2 

ou j represente le nombre imaginaire pur verifiant la relation j 2 = — 1. Posons : 




La solution generale s’ecrit alors : 

y(t) = e~ 5t [Ai e jUDt + A 2 e~ jUDt } 

Sachant que : 



e ju D t _ C os(uj D t) + j sin(wot) 

e -jui D t _ cos (jjj D t) _ j sin(wot) 
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y(t) s’ecrit : 

y(t) = e~ dt [(Ax + A 2 ) cos (u> D t) + j {Ax - A 2 ) sin {w D t)] 
y(t) etant reelle, les nombres complexes Ax et A 2 doivent verifier les relations : 

Ax + A 2 : reel 

Ax — A 2 : imaginaire 

A\ et A -2 sont clone complexes conjugues et peuvent se mettre sous la forme: 



Ax = A' e jv 
A 2 = A' e-i* 



e j(ui D t+ip) e ~j(ui D t+<p) 



y{t) s’ecrit alors : 

y{t) = A'e- st 
soit finalement, en posant A = 2 A' : 

y(t ) = A e~ St cos {wot + g>) 

y{t) peut etre interpretee comrne une fonction periodiciue de pulsation oj u , de phase initiale ip et d’amplitude 
Ae~ St decroissant exponentiellement au cours du temps. On peut definir une pseudo-periode 

T d = — 
lod 

A et ip sont deux constantes d’integration definies a partir des conditions initiales : 



Vo = A cos {ip) 

yo = — A cos{ip) [5 + ujd tan(y>)] 



D’ou l’on tire : 



<P 

A 



— arctan 



y o 

yo + Syo 



Viupyo) 2 + {yo + 8y 0 )i 



Up 



La figure ci-dessous represente les variations de y(t) dans le cas particulier de conditions initiales :yo ^ 
0; y 0 = 0. 
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Remarque : Cas particulier ou <5 = 0: 
L’equation differentielle s’ecrit : 



y + u 2 0 y = 0 

La solution s’exprime dans ce cas : 



y(t) = A cos(wo t + <j>) 

Cette solution est appelee solution harmonique car y(t) est une fonction sinusoi'dale du temps, de pulsation 
wo , de periode To = 27r/wo et dont l’amplitude A et la phase initiale </> sont determinees par les conditions 
initiales : 



d’ou l’on tire : 



Vo = A cos(^) 
y 0 = -w 0 A sin(</>) 



0 

A 



— arctan 



tup y o 
. yo . 





2 



La figure suivante represente les variations au cours du 
1 et y 0 = 0. 



temps de y(t) dans le cas particulier ou : 



y o = 
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temps 

Oscillations non amorties 



3 Equation avec second membre 

3.1 Solution generale 

Soit y(t) la solution generale de Pequation differentielle avec second membre: 

y + 2 5 y + Wq y = A(t) 

Soit yn{t) la- solution de l’equation homogene et soit yp{t) une solution particuliere de Pequation avec second 
membre; yp{t) et yp{t) sont les solutions respectives des deux equations differentielles suivantes : 

Vh + 2 S y H + Wq yH = 0 
jj P + 2 5 y P + Wg y P = A(t) 

Les operations de derivation qui interviennent etant des operations lineaires, l’addition membre membre des 
deux equations differentielles precedentes donne : 

(vh + Vp) + 2 S(y H + y P ) + Wq (y H + yp ) = A(t) 

Ainsi la solution generale peut etre obtenue en faisant la somme de la solution homogene et d’une solution 
particuliere : 



y(t) = yH{t) + yp{t) 

3.2 Cas particulier ou A(t) est constante 

Considerons le cas particulier important ou A(t) = A 0 , Aq etant une constante reelle. L’equation differentielle 
s’ecrit alors : 



On peut verifier que 



y + 25y + ujly = A 0 



An 



y = 
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constitue une solution particuliere de l’equation avec second membre. Selon le cas la solution generale de 
l’equation avec second membre sera : 

• Regime fortement amorti (8 > wq) 



y(t) = A x e s 14 + A 2 e S2t + ^ 

lo 0 

• Regime critique (8 = w 0 ) 

y(t) = (Ai + A 2 t) e St H j 

• Regime pseudo-periodique (8 < ojo) 

y(t) = A e~ St cos(cu£)i + ip) H ^ 

uj 0 

Dans les trois cas consideres les constantes d’integration A± , A 2 , A et cj) son t determinees a partir des 
conditions initiales. Considerons a titre d’exemple le cas particulier y 0 = 0 et yo = 0 . La resolution des 
systemes d’equations obtenus permet de calculer les constantes d’integration et de tracer les variations de y(t) 
pour les trois regimes. Les variations correspondantes de y(t) sont representees par la figure ci-dessous. 




3.3 Cas particulier ou A(t) = A 0 cos (Qt) : 

Nous pouvons verifier que yp(t) = Y 0 cos(f U + 9) constitue une solution particuliere de l’equation differentielle 
avec second membre a condition que l’amplitude Yo et la phase 9 verifient la relation : 

[cap — tl 2 ] Y 0 cos(f2 1 + 9) — 2 SO, sin (fit + 9) = A 0 cos (fit) 

Le developpement des termes en cosinus et en sinus permet d’obtenir : 
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A) 

yj[<4 - n 2 } 2 + 45 2 n 2 



— arctan 



2 SSI 

OJq — fl 2 



La solution complete s’ecrit alors suivant le cas : 
• Regime fortement amorti (5 > wq) 



y(t) = A\e Slt + A 2 e S2t + Y 0 cos(f2 1 + 9) 



• Regime critique (5 = w 0 ) 



y(t) = {A\ + A^t) e St + Y 0 cos (fit + 9) 
• Regime pseudo-periodique (<5 < wq) 



y(t) = A e St cos (coot + <p) + Y 0 cos(f U + 9) 

Comme dans le cas precedent les constantes d’integration A\ , A 2 , A et 9 sont determinees a partir des 
conditions initiales. La figure suivante represente le resultat obtenu dans le cas particulier ou:( yo = 0 , yo = 0 ). 



Regime 

jtransitoirej 

◄ 



Regime permanent 




temps 

Cas ou A(t) est une fonction sinusoi'dale 



On remarque que dans tous les cas la solution homogene tend zero lorsque t augmente, et la solution generate 
s’identifie alors avec la solution particuliere: 

y(t) ~ y P [t) 

Pour cette raison la solution de l’equation homogene est appelee solution transitoire tandis que la solution 
particuliere est appelee solution permanente. 

3.4 Cas ou A(t) est une fonction periodique du temps 

3.4.1 Introduction 

Nous avons etudie dans le paragraphe precedent la solution de l’equation differenticlle dans le cas harmonique, 
c’est-a-dire lorsque le second membre est une fonction sinusoi'dale de la variable t. En considerant le cas de 
fonctions periodiques, nous procederons a une generalisation du cas harmonique. 
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3.4.2 Developpement en serie de Fourier d’une fonction periodique 

Les functions periodiques Une fonction /(f) est dite periodique, de periode T, si pour tout t, f(t+T ) = /(t) 
ou T est une constante positive. La plus petite valeur non nulle de T est appele la periode de f(t). 

Exemples : 



1. La fonction sin(7rt) reprend les memes valeurs pour f = 2,4,6 puisque sin[7r(t + 2)] = sin[7r(t + 4)] = 
sin[7r(t + 6)] = sin(t). Cependant 2 est la periode de sin(7rf). 

2. La periode de sin(flt) ou de cos(fif) est 2tt/S'1. La periode de sin(nfif) et de cos(nfff) est 2ir/n£l. 

3. Une constante a pour periode n’importe quel nombre positif. 

4. Les figures ci-dessous representent deux exemples de fonctions periodiques non sinusoidales. 




Exemples de fonctions periodiques 



Definition des series de Fourier Soit /(t) une fonction periodique de periode T, c’est-dire telle que 
/(t) = f(t + T), la serie de Fourier ou le developpement de Fourier qui correspond /(f) est definie par : 

«‘>-T 

ou n = 2 tt/T, n est un entier positif, a n 
par les expressions suivantes : 



La valeur ao/2 represente la valeur moyenne 
, la plus faible est appele fondamental tandis que les terrnes de pulsation Q, n = nfl plus elevee sont appeles les 
harmoniques. 

Cas des fonctions paires et impaires 

• Une fonction est dite paire si /(— t) = /(t). Dans le cas de developpement en serie de Fourier d’une fonction 
paire,il n’y aura que les termes en cosinus et parfois une constante qui represente la valeur moyenne. II 
est aise de montrer que : 



+ ^ a n cos(nflf) + b n sin (nflt) 

n— 1 

et b n sont les coefficients de Fourier. Ces coefficients sont definis 

/(f) cos(nf!f)cff 
/(f) sin(nflf)df 

de /(f) sur la periode T. Le terme sinusoidal de pulsation fix = ff 




b 



n 



a n 




/(f) cos 




dt 
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• Une fonction est dite impaire si /(— t) = —/(/). Seuls les termes en sinus peuvent tres presents dans un 
developpement en serie de Fourier d’une fonction impaire. II est egalement aise de montrer que : 






b 



n 



o 

4 

T 




f{t) sin 




dt 



Solution de l’equation differentielle La fonction Aft) etant periodique, de periode T, son developpement 
de Fourier s’ecrit : 

OO 

Aft) = -^ + ^ a n cos(nOf) + b n sin(nQi) 

n= 1 

L’equation differentielle s’ecrit alors : 

OO 

y + 25y + u)Qy = — + ^ a n cos(nOt) + b n sin(nflt) 



La reponse permanente ( ou stationnaire ) qui s’iclentifie avec la solution particuliere, pour t suffisamment 
eleve, peut alors etre calculee pour chacune des composantes de l’excitation : oo/2, a n cos(nQt), b n sm(nflt). 
On obtient alors par superposition : 



y(t) = 



a 0 

2fl 2 0 



n—1 



a n cos(f l n t + 0 n ) + b n sin(f \ n t + 9 n ) 
V /(n2- w g ) 2 + 45 2 fi2 



EXERCICES 

Exercice 1 Calculer et representer graphiquement la solution de V equation differentielle homogene : y + dy = 0 
, pour les conditions initiales suivantes : 

• 2/(0) = 1 et y( 0) = 0 

• y(0) = 0 et y{ 0) = 2 

Exercice 2 Calculer et representer graphiquement les solutions des equations differentielles homogenes suiv- 
antes: 

1. y + 5y + 4y = 0 

2. y + 4y + 4y = 0 

3. y + 0.3 y + 4y = 0 

pour les conditions initiales suivantes : 

• 2/(0) = 1 et y{ 0) = 0 

• 2/(0) = 0 et y( 0) = 2 

• 2/(0) = 1 et y{ 0) = 2 
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Exercice 3 Pour les conditions initiates 

2 /( 0 ) = 1 et 2/(0) = 0 

calculer et donner la representation graphique de la solution generate de chacune des equations differentielles 
inhomogenes suivantes : 

1. y + Ay = 5 

2. y + Ay = 5 cos(3 1) 

3. y + Ay = 5cos(2t) 

4. y + 5y + 4y = 5 

5. y + 5y + 4y = 5 cos(3 2) 

6. y + Ay + 4y = 5 

7. y + 4y + 4y = 5 cos(3t) 

8. ?/ + 0.3y + 4y = 5 

9. y + 0.32/ + 4y = 5 cos(3t) 

Exercice 4 Calculer une solution particuliere des equations differentielles inhomogenes suivantes : 

1. 2/ + 42/ = /(t) 

2. y + 4y + 4y = /(i) 

5. y + 5y + Ay = f(t ) 

?/ + 0.3y + 42/ = f(t) 

oil f(t) est une fonction de periode T, definie par : 

f(t) = —a pour—T/2<t<0 
fit) = a pour §<t<T/2 
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Chapitre 1 

Introduction aux equations de 
Lagrange 



1.1 Equations de Lagrange pour une particule 

1.1.1 Equations de Lagrange 

Considerons le cas particulier d’une particule astreinte a se deplacer, sans frottement, sur 
une courbe plane contenue dans le plan xOy. La courbe sur laquelle est astreinte se deplacer la 
particule de masse to, est le lieu des points dont les coordonnees verifient les relation : 



z = 0 

f{x,y) = 0 



La premiere relation correspond au plan xOy . La seconde relation represente l’equation de 
la trajectoire dans ce plan. Ces deux relations definissent. les equations des liaisons appelees 
souvent- liaisons. Le nombre de degres de liberte est egal au nombre de coordonnees qui repre- 
sentent la position de to ( trois dans le cas general) moins le nombre de liaisons ( deux dans le 
cas present ). La particule possede done un degre de liberte. II faut choisir une variable q pour 
reperer sa position. Cette variable est appelee coordonnee generalisee. II est possible d’exprimer 
le vecteur position f de la particule en fonction de la coordonnee generalisee q pr la relation : 
r = f(q). 

Soit F la resultante de toutes les forces agissant- sur la particule. La relation fondamentale 
de la dynamique s’ecrit : 



cPr dv 

r = m—rr = TO — 

dt 2 dt 



ou v = — est la vitesse de la particule. 

Soit 8W le travail fourni par la force F lors d’un deplacement infinitesimal Sr : 



sw = F -Sr 
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Le deplacement Sr peut s’ecrire en fonct-ion de la variation Sq de la coordonnee generalisee 



Q : 

dr 

Sf= —Sq 
dq 

Dans ce cas le travail SW peut se mettre la forme : 

-> dr 

SW = F ■ —Sq 
dq 

On appelle force generalisee conjuguee de q. ou g-composante de la force, la quantite F q 
definie par : 

SW p> dr 
q = ~S^ = "dq 

Par consequent SW s’ecrit : 

SW = F q Sq 

En tenant compte de la relation fondamentale de la dynamique, cette expression peut ega- 
lement s’ecrire : 



dv dr 

SW = m— ■ - —Sq 
dt dq 



D’ autre part 



Sachant que 



d 


_ * dr 


dv dr ^ d 


dr 


dt 


V ' dq. 


dt dq dt 


dq 



d 


dr 


d 


dr 


dt 


dq 


dq 


dt 



dv 

dq 



on obtient 

dv dr d 

dt dq dt 

Le vecteur vitesse v. peut aussi s’ecrire : 



v ■ 



dr 

dq 



— v ■ 



dv 

dq 



D’ou la relation : 



et 



^ dr dr dq dr . 

dt dq dt dq ^ 



dr dv 
dq dq 



dv dr d 

dt dq dt 



v ■ 



dv 

dq 



— v ■ 



dv 

dq 
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1.1 Equations de Lagrange pour une particule 



Sachant que 



et que 



on obtient 



d \ 1 o'! d n 1 _ dv 

~d~q 2 V = Wq 2 V ' V = V ' dj 



d f 1 2 ~| d r 1 J _ dv 
d q 2 V = ~dq 2 V ' V = V ' dq 



civ dr cl f d f 1 



clt dq clt | dq |_2 ? 
L’expression du travail SW peut alors s’ecrire : 



d fl 



dq 2° 



f cl \ d \ 1 2 1 1 d \ 1 2 1 1 

{W=m s b ? Sq 



Si on note T = \mv 1 l’energie cinetique de la masse m , on obtient finalement : 

xu .jd\dT] dT\ 

6 { clt [ dq \ dq J q 

On obtient finalement les deux expressions equivalentes du travail 8W 

( d tar] ari 

\dt[di\--a q r 9 = FA 

On en deduit- l’equation de Lagrange pour un systeme a un degre de liberte : 

cl \dT]dT_ 
clt dq dq q 



1.1.2 Cas des syst ernes conservatifs 

Dans les systemes conservatifs, la force appliquee au systeme derive d’un potentiel U et elle 
s’ecrit : 

F = JJL 

q dq 

L’equation de Lagrange devient alors : 

d rdTl dT dU 



clt I dq I dq dq 



Generalement l’energie potentielle U ne depend pas de la vitesse, c’est-dire que — — = 0. 

dq 

L’equation de Lagrange peut alors s’ecrire : 

d_ \ d(T-U) 1 _ d(T -U) _ 
clt dq dq 
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1.1 Equations de Lagrange pour une particule 
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On introduit la fonction de Lagrange ( ou lagrangien du systeme ) qui est la difference de 
l’energie cinetique et de l’energie potentielle : 

C = T — U 

D’ou la forme de l’equation de Lagrange dans le cas d’un systeme conservatif : 

d_ raci _ dc _ 

dt dq dq 



1.1.3 Cas des forces de frottement dependant de la vitesse 

Equation de Lagrange 

Considerons une situation physique dans laquelle la particule est soumise a des forces de 
frottement de viscosite dont la resultante / est de la forme : 

f = -av 



Pour calculer la force generalisee f q correspondante, nous utilisons la definition du para- 
graphe precedent : 



f q 




= —a 



dr 2 dq 
dq dt 



Cette derniere expression pent se mettre sous la forme : 




d \dC 1 _ dC 
dt dq dq 
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1.1 Equations de Lagrange pour une particule 
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Fonction dissipation 

Calculons le travail SWf fourni par la force de frottement pendant un intervalle de temps 
St pour un deplacement Sr : 

SWf = f ■ Sr = — a v 2 St 

La quantite de chaleur SC) gagnee par le systeme en interaction avec la particule, est telle 
que : 



SQ = a v 2 St 



SQ 



Soit P d = — — la puissance dissipee par les forces de frottement sous forme de chaleur : 
St 

P d = a v 2 

Cette puissance dissipee peut etre exprimee en fonction de q. par : 

2 



Pd = at 



dr 


2 


dr dq 


dt 


= Oi 


dq dt 



= (3 q 



•2 



Par definition, la fonction dissipation est egale a la demi-puissance dissipee : 
La g-composant.e f q de la force de frottement peut alors s’ecrire : 



u = - 



dv 

dq 



L’equation de Lagrange s’ecrit alors : 

d \dC 
dt dq 



dC dv_ 

dq ^ dq 



1.1.4 Cas d’une force exterieure dependant du temps 

Considerons le cas plus general d’un force exterieure dependant du temps agissant sur un 
systeme qui est le siege de forces de frottement qui derivent d’une fonction dissipation V. Soit 
F eq la g-composante de la force exterieure. Dans ce cas l’equation de Lagrange peut s’ecrire 
sous l’une des deux formes equivalentes suivantes : 



d 


~dc 


dc 


dt 


dq 


dq 


d 


' dC 


dC 


dt 


dq 


dq 



Feq~ P q 



dv 

dq 



P 



e,q 
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1.2 Systeme a plusieurs degres de liberte 
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1.2 Systeme a plusieurs degres de liberte 



Dans le cas general d’un systeme a plusieurs degres de liberte, il y a autant d’equations 
de Lagrange que de degres de liberte. Ainsi, si le systeme possede N degres de liberte, il est 
necessaire d’avoir N coordonnees generalisees q, (i = 1, 2, N ) ; nous aurons ainsi N equations 
de Lagrange : 

d 

dt 



dC dC dV ^ .. 1 n 



La q,; — compos ante de la force generalisee exterieure est definie par : 



F = 

1 e,q% 



5W 



Sqi 



Sqi^O 



Dans cette expression 5W represente le travail des forces exterieures resultant d’une variation 
Sqi de la coordonnee q-, telle que les coordonnees q^, soient constantes (Sqj^i = 0). 
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Chapitre 2 

Oscillations libres des systemes a un 
degre de liberte 



2.1 Oscillations non amort ies 

2.1.1 Oscillateur lineaire 

Un systeme oscillant a un degre de liberte est habituellement repere a l’aide d’une coor- 
donnee generalisee q qui est l’ecart par rapport a la position d’equilibre stable. Le mouvement 
vibratoire est dit- lineaire s’il est regi par une equation differentielle harmonique de la forme : 

q + UJ 2 0 q = 0 

Cette equation est appelee equation differentielle de l’oscillateur harmonique simple. 



2.1.2 Energie Cinetique 

Dans le cas d’un systeme a un degre de liberte, constitue d’une masse m dont la position 
est reperee par la coordonnee generalisee q, l’energie cinetique s’ecrit : 



1 2 1 

T = -mv = -m 
2 2 



dr 

dt 



1 

- -m 



dr dq 
dq dt 



1 

~ 2 m 



dr 

dq 



L’energie cinetique d’un systeme a un degre de liberte est fonction de q et q . Elle peut 
s’ecrire sous la forme : 



T = \a(q) q 2 

ou a(q) est une fonction de la coordonnee generalisee q. definie dans le cas etudie par : 



a(q) = m 
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2.1 Oscillations non amorties 
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En faisant un developpement limite de a(q) au second ordre en q. au voisinage de q = 0 , on 
obtient : 

.. 1 ... da 1 d 2 a 2 , 2 

T(q ' q)= 2[ a{0)+ T U J + 2W\J + -\ q 

En limitant P approximation au second ordre, on obtient : 



rri • 2 

T = -ao q 

ou a 0 est une constants egale a a (0) . 

2.1.3 Energie potentielle 

Les oscillations se font autour de la position d’equilibre stable q = 0 caracterisee par 

dU 

——(q = 0) = 0. II est toujours possible , lorsque les ecarts par rapport a la position d’equilibre 
dq 

sont faibles, de faire un developpement en serie de Taylor de U(q ) au voisinage de la position 
d’equilibre q = 0. En negligeant les puissances de q d’ordre superieur a deux, on obtient : 



U(q) = U( 0) 



dU 1 d 2 u 2 

„J + 2 W „J 



q = 0 correspond a un minimum de U (q) pour lequel 

dU d 2 U 

— = 0 et —~y > 0 

d( l q =o d( l q =o 

Si on choisit, l’origine de l’energie potentielle a cett-e position d’equilibre (17(0) = 0) , l’energie 
potentielle U (q) peut s’ecrire sous une forme quadratique : 

U(q) ^ ^ 2 



avec : b n = 



dq 2 „=n 



2.1.4 Equation different ielle 



L’equation de Lagrange s’ecrit : 



d \dC j _ dC _ 
dt [dq \ dq 



Ce qui permet dobtenir l’equat-ion differentielle de l’oscillateur harmonique simple ou : 



d 2 U 

dq 2 n —r 



uj 0 = — = 
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Les oscillations d’un systeme vibratoire s’effectuent autour d’une position d’equilibre stable. 
Pour des oscillations de faible amplitude autour de la position d’equilibre, t-ous les mouvements 
vibratoires peuvent et-re assimiles a des vibrations lineaires et Penergie potentielle peut alors 
etre approximee par une forme quadrat-ique de la coordonnee q. tandis que Penergie cinetique 
peut etre approximee par une forme quadrat-ique en q. 

2.1.5 Resolution de l’equation differentielle de l’oscillateur harmo- 
nique simple 

L’equation differentielle de l’oscillateur harmonique simple s’ecrit : 

q+u 2 0 q = 0 

La solution d’une telle equation est une fonct-ion sinusoi'dale du temps 

q(t ) = A cos (cu 0 1 + (p) 

ou A represent-e l’amplit-ude des oscillations, ip est la phase init-iale. 

II est important de remarquer que la pulsation propre ujo no depend que des elements qui 
constituent le systeme physique et-udie (masse, ressort-, etc...) tandis que l’amplitude A et la 
phase init-iale </? sont calculees a partir des conditions init-iales : 

f = 0) = q 0 

1 q(t = o) = go 

Enfin l’amplitude des oscillations d’un oscillateur harmonique libre ne depend pas du temps. 
De t-elles oscillations sont dites non amorties. 

II faut- neanmoins remarquer qu’au dela d’une cert-aine amplitude la vibration devient non 
lineaire. II s’ensuit d’abord une modification de la periode des oscillations et- ensuit-e un chan- 
gement- de la nature du mouvement. 

2.2 Oscillations libres des systemes amortis a un degre 
de liberte 

Dans le paragraphe precedent, nous n’avons pas tenu compt-e de cert-aines realites physiques. 
En effet-, nous n’avons pas pris en compt-e les forces de frot-t-ement qui sont a l’origine de la pert-e 
d’energie mecanique du systeme sous forme de chaleur. Dans ce paragraphe, nous allons t-enir 
compt-e de ces realites en nous limit-ant- t-out-efois au cas simple ou les pert-es sont dues a des 
frot-t-ements visqueux pour lesquels les forces de frot-t-ement-, qui s’opposent- au mouvement-, sont 
proport-ionnelles a la vitesse . 
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14 



2.2.1 Equation de Lagrange pour les systemes dissipatifs 

Rappelons l’equat-ion de Lagrange associee a un systeme a un degre de liberte dont revo- 
lution au cours du temps se ramene a l’etude de la coordonnee generalisee q 

d_ 

clt 

Fq represente la composante suivant- q de la result-ante des forces generalisees qui ne derivent 
pas d’un potentiel. 

Nous nous int-eressons au cas particulier des forces de frottement definies par la force gene- 
ralisee 



dC _ dC _ 
dq dq q 



F q = ,f q = -fid 

ou 3 est une const-ant-e reelle positive. 

L’equation de Lagrange s’ecrit- alors dans ce cas : 

d dC 
dt dq 




2.2.2 Cas particulier des oscillations de faible amplitude 

Nous avons montre dans le chapit-re precedent que dans ce cas, la fonct-ion de Lagrange 
s’ecrit sous la forme : 

r 1 -2 2 

£ = 2 aq ~ 2 bq 

L’equation different-idle du mouvement s’ecrit- alors : 

a q + bq = ~3q 



C’est une equation different-ielle du second ordre a coefficients constants qui peut se mettre 
sous la forme : 

q + 28q + ulq = 0 

ou 5 est un coefficient positif, appele fact-eur (ou coefficient-) d’amortissement- et- defini par : 



cuo est- la pulsation propre definie par 



s =7r~ 
2 a o 
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2.2.3 Resolution de l’equation different ielle 

La solution de l’equation differentielle depend de la valeur de S par rapport a cu 0 : 

- Si S > ujq . on dit que le systeme est suramorti ou aperiodique. 

- Si 5 = ujq , on dit que Ton a un amortissement critique. 

- Si 5 < ujo , on dit que le systeme est sous-amorti ou pseudoperiodique. 

Cas ou le systeme est suramorti (S > coo) 

La solution de l’equation differentielle s’ecrit dans ce cas : 

q(t) = A 1 el~ S -^ Z ^ t + A 2 el- 5+ ^^\ t 

Ai et A 2 sont des constantes d’integration definies par les conditions initiates. La figure 
ci-dessous represente q en fonction du temps dans le cas particulier ou q( 0) = go et g(0) = 0. 
q(t) est une fonction qui tend exponentiellenrent (sans oscillation) vers zero. 




Regime fortement amorti : variation de q en fonction du temps 



Cas de l’amortissement critique (5 = cu 0 ) 

La solution generate de l’equation differentielle est de la forme : 

q(t) = (A\ + A 2 1) e 5t 

Dans le cas particulier ou g(0) = go et g(0) = 0, 

q(t) = g 0 (1 + St) e~ 5t 

q{t ) est encore une fonction qui tend vers zero sans oscillation lorsque le temps augmente. 
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Amortissement critique : variation de q en fonction du temps 



Cas ou le systeme est sous-amorti (5 < u; 0 ) 

La solution generale de l’equation differentielle est de la forme : 

q[t) = A e~ St cos (uja t + 4>) 



avec uja = a/cUq — 8 2 ; A et 0 sont deux constantes d’integration determinees a partir des 
conditions initiales. Dans le cas particulier ou g(0) = go et g(0) = 0, on obtient : 



A 



UJ 0 

1 

UJ A 



-Qo 



arctan — 




Systeme faiblement amorti : variation de q en fonction du temps 
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2.2 Oscillations libres des systemes amortis a un degre de liberte 



IT 



La figure represente les variations de q(t) au cours du temps. On remarque que q(t) est 
enveloppee par les deux fonctions exponentielles : q 0 e~ St . 

Le lieu des maxima est obtenu en resolvant q(t) = 0. Ce qui donne : 

£ 

tanf/u^t + </>) = = tan (</>) 

UJa 

d’ou Ton tire l’instant t n correspondant au n-ieme maximum : 

2vr 

t n — Tl 

UJA 



Les maxima de q(t) sont separes par des intervalles reguliers egaux a 




UJ A 



Ta est appelee la pseudo-periode. On remarque que , en plus de la diminution de F amplitude 
des oscillations au cours du temps, Fun des effets de Famortissement est Faugmentation de la 
periode des oscillations. Pour des systemes faiblement amortis (5 « uj 0 ), on peut remarquer 

27T 

que uja — <^o et que la pseudo periode est peu differente de la periode propre : Ta — T 0 = — . 

cu 0 



2.2.4 Exemples 

Systeme mecanique en translation 

Amortisseur mecanique Un amortisseur mecanique est constitue d’un element mobile a 
l’interieur d’un recipient contenant un fluide visqueux. 



a 



















Amortisseur 



La force de frottement / agissant sur la partie mobile reperee par xq. est donnee par 

fx = -a (;r i - ;r 2 ) 

ou (.iq — i' 2 ) represente la vitesse relative des deux elements qui constituent Famortisseur. 
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Equation differentielle du mouvement Considerons le cas d’une masse m oscillant ver- 
ticalement et reliee a un bati fixe par un ressort de raideur k et un amortisseur de coefficient 
de frottement visqueux a. Reperons par x l’ecart de la masse m par rapport a la position 
d’equilibre.L’equation de Lagrange s’ecrit- : 

d_ lac 

clt dq 

Sachant que T = \mx 2 et que U = \kx 2 , la fonction de Lagrange se met sous la forme 

r 1 • 2 1 7 2 

L = -mx — -Kx 
2 2 

On obtient l’equation differentielle du mouvement 

mx + ax + kx = 0 

Cette equation pent s’ecrire sous la forme d’une equation differentielle du second ordre a 
coefficients constants 



dC 

dq 



= —ax 



x + 2 5 x + (jj^x = 0 

ou cuq et 5 sont des const-antes positives appelees respectivement la pulsation propre et le 
facteur d’amortissement- ; elles sont donnees par les relations suivantes 



5 



a 

2 m 



et cuq 




Systeme mecanique en rotation 

Considerons un pendule simple const-itue par une masse ponctuelle m reliee a un axe de 
rotation fixe par une tige de masse negligeable, qui effectue des oscillations de faible amplitude 
dans un plan vertical. Ce disposit-if est un systeme a un degre de liberte dont on repere la 
position par Tangle 0 par rapport a la vert-icale. Le systeme est immerge dans un fluide dont les 
forces de viscosite se ramenent a une force s’exergant sur la masse m et donnee par la relation 

/ = ~av 



L’equation de Lagrange qui regit le mouvement d’un tel disposit-if s’ecrit : 

cl 
clt 

L’equat-ion differentielle du mouvement s’ecrit- alors : 

ml 2 9 + a l 2 9 + m g 1 9 = 0 



DC 

dq 



dC 

dq 



= —al 2 6 
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9 + 259 + co 2 0 9 = 0 

ou la pulsation propre cu 0 et le facteur d’amortissement 5 sont respectivement donnes par : 



UJ 0 




et 



5 



a 

2 m 
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Chapitre 3 

Oscillations forcees des systemes a un 
degre de liberte 



3.1 Equation different ielle 

Rappelons la forme generale de l’equat.ion de Lagrange pour les systemes a un degre de 
liberte : 



d \0C 
dt dq 



dC dV_ 

, , r „ . -C next 

oq oq 



ou F qext est la force generalisee associee a F ext et ou la fonction dissipation est V = 

2^ 2 ' 

Pour les oscillations de faible amplitude, la fonction de Lagrange pouvait se mettre sous une 
forme quadrat-ique de q et q 

r 1 . 2 ! ? 2 

£ = 2 °° V _ 2° q 

D'ou l’equation differentielle du mouvement 



°o Q + P Q + b 0 q — F qext 

Cette equation peut, se mettre sous la forme d’une equation differentielle du second ordre a 
coefficients constants, avec second membre 

q + 28q + ulq = A(t ) 



avec 



5 



A 

2ao ’ 



OJ 0 




et A(t) 



Fqext 

do 
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3.2 Systeme masse-ressort-amortisseur 




m 



x 



F(t) 



Systeme masse-ressort-amortisseur 



Considerons l’exemple mecanique de la figure ci-dessus soumis a une force exterieure F(t) 
appliquee a la masse m. Calculons la force generalisee F x conjuguee de la coordonnee x. Pour 
cela nous pouvons utiliser l’une des deux methodes suivantes : 

- Soit calculer le travail dW de la force F(t) pour une variation dr de son point d’application 

dW = F ■ dr = F dx 



On en deduit la .r-composante de la force exterieure 

^ = FW 

- Soit utiliser la definition de la force generalisee 

— * f)r 

F . = F -is = m 

L ’equation differentielle du mouvement s’ecrit alors 



avec : 



x + 2 6x + (jj^x = A(t) 



5 = a/2m , ujq 




A(t ) = F(t)/m 
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3.3 Solution de l’equation differentielle 

La solution de cette equation differentielle du second ordre est egale a la somme de la solution 
de 1’ equation sans second membre (ou solution homogene) x pit) et d’une solution particuliere 
de l’equation avec second membre x pit) : 

x(t) = x H (t) + Xp{t ) 

Nous avons deja etudie l’equat-ion sans second membre Xpit) et nous savons que cette 
solution contient dans t.ous les cas le terme exponentiel e~ 5t . Apres un temps t superieur a 3/(5 ou 
4/5, le terme e 5t devient tres petit et la solution homogene est alors pratiquement nulle. II ne 
subsistera que la solution particuliere de l’equation avec second membre. L’intervalle de temps 
pendant lequel la solution homogene est non negligeable est appele le regime transitoire. A la 
fin de ce regime transitoire commence l’intervalle de temps pour lequel la solution homogene 
est quasi- nulle et pour lequel la solution x(t) ~ x p (t) ; ce regime est appele regime permanent 
ou stationnaire. 



3.3.1 Cas particulier ou A(t) = Aq cos(f2t) 

a) Calcul de la solution permanente l’aide de la methode des nombres 
complexes 

Pour t suffisamment grand, nous pouvons considerer que la solution transitoire s’est- annulee 
et que la solution x(t) s’identifie alors avec la solution particuliere : x(t) ~ Xp(t). Par commodite 
de notation l’indice p est sous-entendu dans ce qui suit. La methode des nombres complexes 
permet de calculer aisement la solution stationnaire. 

Soit- le deplacement complexe represente par le nombre complexe X = avec X = 

X 0 e J:p . Nous pouvons considerer, en outre, que A(t) = A 0 cos (Qt) constitue la partie reelle du 
nombre complexe A = Aq e j ' H . L’equation differentielle se transforme en une simple equation 
algebrique en fonction de l’amplitude complexe A : 

[(wg -fl 2 ) + j25Q] X = A 0 

dont- la solution est : 



— (uj 2 0 -n 2 )+j2sn 

D’ou Ton tire l’amplitude X 0 et la phase <p : 



y / (a;2-p2)2 + 4(5 2 fi2 

26fl 

if = — arctan — » — 

cUq — f V 
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b) Etude des variations de l’amplitude et de la phase en fonction de 
la pulsation de l’excitation 

dX 

Le maximum de l’amplitude est obtenu pour la valeur de Q qui annule — — . 

ail 

II existe un maximum a la pulsation Q R = a/cUq — 2 8 2 seulement si l’amortissement est 
suffisamment faible pour que 8 < \[2- A cette pulsation appelee pulsation de resonance, on 

dit que le systeme entre en resonance et l’amplitude X 0 est maximale ; elle vaut : 



X, 



Ac 



0 max 



2 Sy/ul - 8 2 

La figure representant les variations de X 0 en fonction de la pulsation d’excitation 12 est 
appelee courbe de resonance en amplitude. On remarque qu’a la pulsation uj 0 , le dephasage ip 



7T 



est egal a — — , et qu’a la resonance ip = — arctan 

At 



\Ju%- 2 8‘ 2 
8 




c) Etude de la resonance pour les faibles amortissements 

Dans le cas des faibles amortissements ( 8 « oj 0 ), la frequence de resonance est tres peu 
differente de la pulsation propre, Q R ~ cu 0 . Dans ce cas, l’amplitude de vibration a la resonance 
Aomax est egale a : 

X - 
av 0 max n e 

Zouj o 

Pour les faibles amortissements, X (} rnax est done inversement proportionnel a 8. 



d) Etude de la vitesse 

En notation complexe, la vitesse s’ecrit : 

V(t) = ^ = jVLX = 
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ou 1’ amplitude complexe de la vitesse est definie par 



X = jVtX 



jQA 0 

(u 2 0 -W)+j25tt 



L’etude des variations de l’amplitude de la vitesse en fonction de la pulsation d’excitation 
montre que, quelle que soit la valeur de <5, la resonance en vitesse est obtenue pour Q = u (J 
(voir figure ci-dessous). La valeur maximale de l’amplitude de la vitesse vaut dans ce cas : 



di-max — -A- (n^o ) 



Aq 

25 




rc/2 



0 



-n/2 




Courbe de resonance de la vitesse Dephasage ip de la vitesse en fonction de Q 



e) Bilan energetique 

Soit Pp(t) la puissance instantanee fournie par la force exterieure F(t) au systeme. En 
regime permanent, on obtient : 

Pp(t) = F(t) x(t ) = F 0 X 0 cos (fit) cos (fit + ip) 

Soit < Pp > la valeur moyenne sur une periode de Pi-(t) : 

< Pf >= t^FqXq cos (ip) 

En tenant conrpte de Pexpression de X 0 en fonction de F 0 , on obtient : 

1 • 9 

< Pf >= 2 aX o 

Comparons cette valeur a la valeur moyenne < Pp > de la puissance dissipee par les forces 
de frot.tement de viscosite. La valeur instantanee de cette puissance dissipee s’ecrit : 

Po(t) = ax 2 = aX 2 cos 2 (fit + ip) 

D’ou l’on tire la valeur moyenne sur une periode : 



http:/ / dj elouah.ifrance.com 



3.3 Solution de l’equation differentielle 



25 



1 ' 9 

< Pd >= 2 aX o 

L’etude des variations de la valeur moyenne de la puissance < P >=< Pp >=< Pd > en 
fonction de la pulsation d’excitation montre que la valeur maximale de la puissance moyenne 
est obtenue pour = co 0 quelle que soit la valeur de S. La valeur maximale de la puissance 
moyenne dissipee ou fournie vaut dans ce cas 

F 2 

<P> = — 

\ i /max 0 

2a 

La figure ci-dessous represente les variations, en fonction de Q. de la puissance moyenne 
dissipee par les forces de frottements ( ou de la puissance moyenne fournie par la force exterieure 

)• 




Courbe de resonance pour la puissance 



f) Bande passante 

On definit par bande passante, la bande des pulsations autour de 0 = cu 0 pour lesquelles 
< P >>< P > m ax /2. Les deux pulsations et 0 2 ,situees de part et d’autre de la pulsation 
ujo et pour lesquelles < P >=< P > max /2, sont appelees pulsations de coupure. La bande 
passante B s’ecrit : 



B = — 

Le calcul de B consiste a rechercher les deux pulsations pour lesquelles < P >=< P > max / 2. 
On obtient l’expression de la bande passante B : 
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B = n 2 -n 1 = 2S 

g) Coefficient de qualite d’un oscillateur 

Le coefficient de qualite d’un oscillateur est defini par le rapport de la pulsation propre cu 0 
a la largeur de bande B : 



Q = 



cu 0 

~B 



3.3.2 Cas d’une excitation periodique 



Nous avons etudie dans le paragraphe precedent la reponse d’un systeme vibratoire a une 
excitation sinusoi'dale dite excitation harmonique. En pratique, les excitations mecaniques ne 
sont pas toujours parfaitement sinusoi'dales ; elles sont souvent periodiques. En considerant le 
cas d’excitat.ions periodiques, nous procederons a une generalisation du cas harmonique. 

Soit une excitation periodique appliquee a un systeme amorti a un degre de liberte. L’equa- 
tion differentielle qui regit ce systeme s’ecrit : 

q + 28q + culq = A(t ) 



La fonction A(t) etant periodique, de periode T, son developpement de Fourier s’ecrit : 



m 



a o 
~2 



OO 

+ a n cos (nut) + b n sin(ncut) 

n= 1 



L’equat.ion differentielle s’ecrit alors : 



q + 2 5 q + cug q 



a o 
~2 



OO 

+ a n cos (nut) + b n sin(ncut) 

n = 1 



La reponse permanent^ ( ou stationnaire ) qui s’identifie avec la solution particuliere, pour 
t suffisamment eleve, peut alors etre calculee pour chacune des composantes de l’excitation : 
Oq/2, a n cos(nujt), b n sm(nujt). On obtient alors par superposition : 



q(t) 



a o 



n = 1 



a n cos {u n t A'lpn) + b n sin (u n t + 'ipj 




A8 2 cu 



2 

n 



3.4 Impedance mecanique 

3.4.1 Definition 

Considerons un systeme mecanique soumis a une force sinusoi'dale F (t) = F () cos (At). En 
regime permanent, le point d’application de cette force se deplace avec une vitesse v (t) = 
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Vo cos (Qt + cf)) . On appelle impedance mecanique d’entree du systeme mecanique, le rapport 
des amplitudes complexes de la force F et de la vitesse v 



3.4.2 Impedances mecaniques 

Amort isseur 

Dans le cas d’un amortisseur, la force appliquee est reliee a la vitesse par 

F = av 

On en deduit l’impedance complexe d’un amortisseur 

Z_ a = a 



Masse 

Dans le cas d’une masse, la relation fondamentale de la dynamique s’ecrit 

civ 

r = m— 
clt 

On en deduit l’impedance complexe d’une masse 

Z_ m = jmtt = mQ e 7 2 

Ressort 

Dans le cas d’un ressort de raideur k, la force appliquee / appliquee au ressort s’exprime en 
fonction de l’allongement par 

f = kx 

On en deduit l’impedance complexe d’un ressort 

k .k k _a 2l 

- k = ]n = ~ J n = n e 

3.4.3 Puissance 

La valeur moyenne, sur une periode, de la puissance fournie est 



< Pf >= 7 >F q X 0 cos (0) = ^ Re (Z E ) X { 



http:/ / dj elouah.ifrance.com 






3.4 Impedance mecanique 



28 



3.4.4 Applications 

Systeme mecanique resonant 

Soit mi systeme mecanique constitue d’un ressort de raideur k, d’un amortisseur de coef- 
ficient de frottement visqueux a et d’une masse m soumise a une force sinusoi'dale F (t) = 
F 0 cos (fit). L’impedance d’entree de ce systeme est 



Z E = a+j 




k 

n 



A la resonance = cuq = 



, le module de l’impedance est Z E = a. Lorsque la pulsa- 



tion — > oo, l’impedance Z_ E — jrnQ. 





Module de l’impedance d’entree Amplitude de la vitesse 

Systeme antiresonant 

Considerons un circuit constitue par un ressort de raideur k dont une extremite est reliee 
a une masse rn et dont l’autre est soumise a une force sinusoi'dale F (t). Soit x le deplacement 
de la masse m et soit y le deplacement du point d’application de la force F (t). Pour calculer 
l’impedance d’entree de ce systeme, nous devons d’abord ecrire les equations differentielles du 
mouvement : 



mx = k (x — y) 

F = k (x — y) 

En utilisant la notation complexe, on obtient- l’impedance d’entree : 
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& = f = -» 



km 






La pulsation d’antiresonance est cuo = \ — • Lorsque Q = cuo, la vitesse Y est nulle t.andis 

V m 

que le module de l’impedance est oo. Lorsque la pulsation — > oc, l’impedance Ze — > 0. 





Module de l’impedance d’entree Amplitude de la vitesse 
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Chapitre 4 

Oscillations libres des systemes a deux 
degres de liberte 

4.1 Introduction 

Les systemes qui necessitent deux coordonnees independantes pour specifier leurs positions 
sont appeles systemes a deux degres de liberte. 

Exemples 




Figure 1 



Figure 2 



Figure 3 



- Figure 1 Si les masses mi et m 2 sont astreintes a se deplacer verticalement, 2 coordonnees 
x 1 et a : 2 sont necessaires pour specifier la position de chaque masse a chaque instant. 

- Figure 2 Si la masse M est astreinte a se deplacer dans un plan vertical, deux coordonnees 
sont necessaires pour specifier la configuration du systeme. L’une de ces coordonnees 
peut etre le deplacement x qui correspond a la translation verticale de la masse. L’autre 
coordonnee peut etre le deplacement angulaire 9 pour t-enir compte de la rotation de la 
masse. Ces deux coordonnees sont independantes l’une de l’autre. 

- Figure 3 Dans le cas du double pendule, deux coordonnees sont necessaires pour specifier- 
la position des masses m ] et to 2 . Plusieurs choix sont pourtant possibles, en effet on peut 
choisir (x u x 2 ) ou (y u y 2 ) ou (9 1 ,9 2 ). 
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II est possible de specifier la configuration d’un systeme a l’aide de plusieurs ensembles de 
coordonnees independantes ; un ensemble quelconque de ces coordonnees est appele coordonnees 
generalisees. II y a autant d’equations de Lagrange que de degres de liberte ou de coordonnees 
generalisees. Pour P etude des systemes a deux degres de liberte, il est necessaire d’ecrire deux 
equations differentielles du mouvement que l’on peut obtenir a partir des equations de Lagrange 



f cl 


'dc 


DC 


1 clt 


dqi 


dqi 


1 d 


'dC 


dC 


{ Jt 


_dq 2 


1 

Qjl 

to | 



4.2 Systemes a deux degres de liberte 

4.2.1 Systeme masses-ressorts en translation 







x 2 



Considerons le systeme ci-dessus, constitue de deux masses m iet m 2 reliees respectivement 
par deux ressorts de raideur k \ (it k 2 a deux batis fixes. Les deux masses sont reliees par un 
ressort de raideur K. Ce ressort est appele ressort de couplage. 

Equations differentielles du mouvement 

Les equations du mouvement pour ce systeme a deux degres de liberte peuvent etre obte- 
nues a partir des equations de Lagrange pour chaque coordonnee X\ (t) et x 2 (t). Soit T et U 
respectivement Penergie cinetique et Penergie potentielle 
T = \ mi x\ + \m 2 x\ 

U = \ki x\ + \K (xi - x 2 ) 2 + \k 2 x% 

U = \ (ki + K ) x\ + l ( k 2 + K ) xl - Kxix 2 
Le lagrangien L = T — U s’ecrit alors 

C = i m i xl + ^ m 2 xl - i (Jfei + K) x\ — ^ (k 2 + K) x\ + Kx rx 2 
Les equation de Lagrange s’ecrivent 
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f d 


'dc 


oc 


I dt 


dxi 


dx\ 


1 d 


'DC 


dc 


[ dt 


Q 3 

to 


dx 2 



D'ou le systeme d’equations different ielles du mouvement 



f m\X i + {h + K ) X\ — Kx 2 = 0 
\ m 2 x 2 + (k 2 + K) x 2 - ay = 0 

Les termes —Kx 2 et —Kx\ qui apparaissent respectivement dans la premiere et la seconde 
equation sont appeles termes de couplage, et les deux equations differentielles sont dites cou- 
plees. 



Resolution des equations differentielles 

Les deux solutions de ces deux equations differentielles sont des fonct.ions periodiques et sont 
composees de deux fonctions harmoniques de pulsations differentes et d’amplitudes differentes. 
Supposons que l’une de ces composantes harmoniques s’ecrive 

xi(t) = Ai cos(ut + 0) 
x 2 (t ) = A 2 cosioot + (f)) 

ou Ai, A 2 et o sont des constantes et to l’une des pulsations propres du systeme. La substi- 
tution de x | et x 2 dans le systeme d’equations differentielles donne 

f [ki + K — mi w 2 ] Ai — K A 2 = 0 
[ — K A\ + [k 2 + K — m 2 a? 2 ] A 2 = 0 

Ce qui constit-ue un systeme d’equations lineaires homogenes dont- les inconnues sont A\ et 
A 2 . Ce systeme admet une solution non identiquement nulle seulement si le determinant A(cu) 
des coefficients de Ai et A 2 est egal a zero. 

[h + K — TYl\ cu 2 ] — K 

- K [k 2 + K - m 2 oj 2 } 

Le determinant A(cu) est appele determinant caracteristique. L’equation A (ou) = 0 est 
appelee l’equation caracteristique ou equation aux pulsations propres. Elle s’ecrit 



A(w) = 



k\ + K — miuj 2 k 2 + K — m 2 oj 2 — K 2 = 0 



ou encore 



4 2 

L 0 ~ LO 



k\ + K , k 2 + K 



mi 



m 2 



k\k 2 + k\K + k 2 K 
m± m 2 



= 0 



Cette equation est une equation quadrat-ique en u qui admet deux solutions reelles positives 
oji et cu 2 appelees les pulsations propres du systeme 
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Cet exemple montre qu’il y a en general deux pulsations propres dans un systeme a deux 
degres de liberte. Chacune des coordonnees, x i ct x 2 , possede deux composantes harmoniques 
de pulsations uqet cu 2 



xi = A u cos(aq t + fa) + A 12 cos(c u 2 t + fa) 
x 2 = A 2 1 cosfuqt + fa) + A 22 cos(c o 2 t + fa) 

ou An, A 12 , A 2 i . A 22 , 0] et fa sont des constantes. Le terme de plus basse frequence 
correspondant- a la pulsation uq est appele le fondamental. L’autre terme, de pulsation u 2 , est 
appele harmonique. 

Les doubles indices sont utilises pour les amplitudes des differentes composantes harmo- 
niques ; le premier indice se refere a la coordonnee et le second a la pulsation. Par exemple A 12 
est l’amplitude de .i'i (t) a la pulsation x 2 - 

Lorsque A i2 = A 22 = 0, aqet x 2 correspondent a la premiere solution particuliere sont des 
fonctions sinusoi'dales, en phase, de pulsation u q ; on dit que le systeme oscille dans le premier 
mode. Dans ce cas 



xi = An cos(aqf + fa ) 
x 2 = A 2 1 cos(cuit + fa) 



Lorsque A n = A 2 \ = 0, aqet x 2 correspondent a la seconde solution particuliere et sont des 
fonctions sinusoi'dales, en opposition de phase, de pulsation ui 2 ; on dit que le systeme oscille 
dans le second mode. Dans ce cas 



X 1 = A 12 cos (cu 2 t + fa) 

x 2 = A 22 cos(u 2 1 + (j) 2 ) 



Etudions les particularites de ces deux solutions particulieres : 
- La premiere solution particuliere s’ecrit : 



Xi = An cos(aqf + 0 1 ) 
x 2 = A 2 1 cos(o;it + fa) 

Xi et x 2 doivent verifier le systeme d’equations differentielles, ce qui donne 

f [h + K - mi fa] A n - K A 2 i = 0 
— A An + [k 2 + A — uq <^i] A 2 i = 0 

Ces deux equations permettent d’obtenir le rapport des amplitudes dans le premier mode 
ou fondamental 

A 2 i ki + K — mi(jj\ K 

^ 1 An K k 2 + K - m 2 uj\ 
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- La seconde solution particuliere s’ecrit : 

xi = A 12 cos(c o 2 t + <j) 2 ) 

X 2 = A 2 2 cos(u 2 t + <j) 2 ) 

x i et x-2 doivent verifier le systeme d’equations differentielles, ce qui donne 

f [k\ + K — nil w 2 ] -4-12 — L’ A 22 = 0 
— K A\ 2 + \k 2 + K — m 2 1 u|] A 22 = 0 

Ces deux equations permettent d’obtenir le rapport des amplitudes dans le second mode 
ou harmonique 

A 2 2 k\ + K — ?7ii uj\ K 

A\2 K k 2 + K - m 2 u\ 

- La solution generate (xi,x 2 ) est une combinaison lineaire de ces deux solutions parti- 
cuieres. aqet x 2 s’ecrivent alors 

x 1 = An cos (u>it + <f>i) T A\2 cos {u 2 t -\- 0 2 ) 
x 2 = /i 1 An cos (uit + 0 i) 4 - /i 2 A12 cos (uj 2 t + (j) 2 ) 

ou An, A 12, 4 >i et (j ) 2 sont des constantes d’integration dont les valeurs sont fixees par les 
conditions initiales. 

4.2.2 Cas particulier de deux oscillateurs identiques 

Calcul des constantes d’integration 

Considerons le cas particulier de deux oscillateurs identiques tels que rrii = m 2 = m et 
A- 1 = A- 2 = k. Dans ce cas les pulsations propres sont respectivement egales a 




Les rapports d’amplitudes correspondant. a ces pulsations sont respectivement // , = +1 et 
Ih = - 1 - 

Soit- x 10, X 20 , xiq et x 2 o les valeurs initiales respectives de xi, x 2 , xi et i ; 2 . Tenant compte de 
ces conditions initiales, on obtient le systeme d’equations suivant qui permet de determiner les 
constantes d’integration An, Ai 2 , 4 >i et (j) 2 

An cos(0 1 ) + Ai 2 cos (<j) 2 ) = aqo 
An cos (</>i) - A 12 cos(</> 2 ) = x 20 
-Ui An sin((j)i) - u> 2 A u sin(0 2 ) = xi 0 
-loi An sin^) +cu 2 4i 2 sin(0 2 ) = x 2 o 

Les solutions de ce systeme d’equations sont 
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ou encore 



x io + %20 
2 cos (00 



et 



x 10 — x 20 
2 cos (0 2 ) 



&10 + x 20 
2 oji sin(0 1 ) 



et 



Ai2 — 



x 20 ~ x 10 
2 cu 2 sin (^2 ) 



1. Considerons le cas particulier suivant x\q = x 2 q = Xq et x w = x 20 = 0 ; on obtient dans 
ce cas cj) 1 = 02 = 0 , A 12 = Oet A n = x 0 ; d’ou 



X\ = x 0 cos(u>it) 
x 2 = Xq cos(u>it) 



Pour ces conditions initiales particulieres, les deux masses oscillent en phase a la meme 
pulsation uj\. On dit que le systeme oscille dans le premier mode. 




2. Considerons un autre cas particulier pour lequel xio = — x 2 q = xq et x \ o = i’ 2 o = 0 . On 
obtient dans ce cas 0i = 0 2 = 0, An = 0 et A 12 = x 0 ; d’ou 

X\ = xq cos(cu 2 t) 
x 2 = — Xq COs(u 2 t ) 



On dit que le systeme oscille dans le second mode car les deux masses oscillent en oppo- 
sition de phase avec le meme pulsation u 2 . 
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3. Considerons enfin le cas particulier suivant x w = x 0 , x 2 o = 0 et x 10 = x 20 
0i = 4*2 = 0; ^ii = -4-12 = xq/2. Les solutions s’ecrivent alors sous la forme 



0 ; d’ou 



x l(t) = cos ( W 1 t) + ^ cos (cu 2 t ) 

X 2 (t) = COS (cui t) — cos (cu 2 t ) 

Les solutions ne sont plus des fonctions purement sinusoi'dales du temps mais des combi- 
naisons lineaires de deux fonctions sinusoi'dales de pulsations respectives uj\ et ui 2 . x± et 
x 2 peuvent s’ecrire sous la forme 



La figure suivante 
(c’est-a -dire si K 



X\ ( t ) = Xq cos 

x 2 ( t ) = xq sin 



UJ 2 - UJ 1 
. 2 
UJ 2 — LUi 



t J COS 
t ) sin 



(X 2 +UJ 1 
v 2 1 
0J 2 + LOi 



represente le resultat obtenu dans le cas ou cui est tres different de co 2 

» k ). 




Si uji est peu different de cu 2 (c’est-a -dire si K « k). on observe un phenomene de bat- 
tement (voir figure ci-dessous). 
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Coordonnees principales 

Considerons les coordonnees p\ et p 2 obtenues a partir des coordonnees x \ et x 2 par les 
relations 



Pi 

P‘2 



xi + X2 
2 

X\ — X2 
2 



Tenant compte des expressions de x i et x 2 et des valeurs particulieres de // , et p 2 I )01ir 
l’exemple etudie, on obtient 



Pi = — cos (uit) 

P2 = y cos (u 2 t) 

On remarque que, quelles que soient les conditions initiales, p\ et p 2 sont des fonctions 
purement sinusoi'dales du temps de pulsations respectives uj i et cu 2 . Ces coordonnees particulieres 
sont appelees coordonnees principales. On pent verifier que le systeme d’equations differentielles 
qui regit le mouvement du systeme considere s’ecrit sous la forme de deux equations decouplees 

pi + u\ pi = 0 
p 2 + ul p 2 = 0 

Les relations inverses suivantes 



xi=pi+ p 2 
x 2 =Pi- p 2 

permettent d’obtenir les coordonnees aqet x 2 a partir des coordonnees principales p\ et p 2 . 
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4.2.3 Pendules couples 

Considerons le cas de deux pendules simples identiques couples par un ressort de raideur K 
et qui effectuent des oscillations de faible amplitude reperees par les angles 9 1 et 0 2 . 




Etablissons tout d’abord les equations differentielles du mouvement dans le cas des oscilla- 
tions de faible amplitude. II est aise de montrer que Penergie cinetique et Penergie potentielle 
s’ecrivent sous les formes quadratiques suivantes 

T = bm/ 2 6 l + ^ml 2 0 2 

U = \ [Kl 2 + mgl ] d\ + § [. Kl 2 + mgl \ Q\ - Kl 2 6 10 2 

On remarque la presence du terme de couplage —Kl 2 0i0 2 dans Pexpression de Penergie 
potentielle. Comme dans Pexemple precedent, on dit que le couplage est elastique. Si le terme 
de couplage n’existe que dans Pexpression de Penergie cinetique, on dit que le couplage est de 
type inertiel. 

Les equations de Lagrange permettent d’obtenir les equations differentielles du mouvement 

m/ 2 0! + [Kl 2 + mgl] 9 1 - Kl 2 9 2 = 0 
—Kl 2 6\ + ml 2 9 2 + [Kl 2 + mgl] 9 2 = 0 

En Pabsence d’amortissement la solution de ce systeme d’equations differentielles est de la 
forme 



6i(t) = Ai cos(c ot + 4>) 

0 2 (t) = A 2 cos(cut + 4>) 

Ces deux expressions doivent satisfaire le systeme d’equations differentielles, d’ou 

[Kl 2 + mgl - ml 2 oj 2 ] A 1 - Kl 2 A 2 = 0 
—Kl 2 Ai + [Kl 2 + mgl - ml 2 oj 2 ] A 2 = 0 

Ce systeme d’equations admet des solutions non nulles seulement si cu est solution de liqua- 
tion aux frequences 
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[. Kl 2 + mgl — ml 2 uj 2 ~\ 2 — K 2 l 4 = 0 
D’ou Ton tire l’expression des pulsations propres u \ et cu 2 




Les solutions du systeme d’equat.ions differentielles sont done 

0 1 = An cos(cuit + (f>i) + Ai 2 cos(uj 2 t + <^2) 

0 2 = -d-21 COS (cUit + 0 X ) + A 22 COS (u 2 t + (j) 2 ) 

Pour calculer les rapports des amplitudes dans les modes, on suppose que le systeme oscille 
soit dans le premier mode soit dans le second mode. Dans le premier mode, on obtient le systeme 

[Kl 2 + mgl — ml 2 (jjf\ — Kl 2 y u 1 = 0 
— Kl 2 + [Kl 2 + mgl — ml 2 uj \ ] /.q = 0 

Dans le second mode, on obtient 

[Kl 2 + mgl — ml 2 uj^\ — Kl 2 n 2 = 0 
— Kl 2 + [Kl 2 + mgl — ml 2 ix> 2 ] /j 2 = 0 

Tenant compte des expressions de uj 1 ct cu 2 on obtient les valeurs du rapport des amplitudes 
dans les modes /i 1 = +1 et /i 2 = —1. Les solutions du systeme d’equations differentielles 
s’ecrivent alors 



9 1 — An cos ( Uit + (f>i ) + Ai 2 cos (oj 2 t + 0 2 ) 
0 2 = -4ll COS ( OJit + (fil) - Ai 2 cos (oj 2 t + (j) 2 ) 
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Chapitre 5 



Oscillations forcees des systemes a 
deux degres de liberte 

5.1 Equations de Lagrange 

Soit- un systeme a deux degres de liberte, soumis a des forces qui derivent d’un potentiel, a 
des forces de frottement de viscosite et des forces exterieures. Si les coordonnees generalisees 
sont q 1 et q 2 , les equations de Lagrange s’ecrivent : 

d_ \dC] _d£ L OV 

clt dqi dqi dqi qi 

d_ \dC] _dC dV 

dt _dq 2 _ dq 2 + dq 2 92 

Dans cette expression F qi et F, n sont les forces generalisees conjuguees des coordonnees 
generalisees respectives q\ et q 2 . Elies sont respectivement definies par 

- F qi = 5qi ^ 0 , dans cette expression 8W\ represente le travail des forces exterieures pour 

< 5^2 =0 

une variation dqi de la coordonnee q\. lorsque 5q 2 = 0. 

- F q2 = 'jdd Sqi=Q , dans cette expression SW 2 represente le travail des forces exterieures pour 

Sq2^=0 

une variation bq 2 de la coordonnee q 2 , lorsque Sq 2 = 0. 

5.2 Systeme masses-ressorts-amortisseurs 

Pour etudier les particularites des oscillations forcees des systemes a deux degres de liberte, 
etudions le systeme symetrique suivant soumis a une force horizontale F. appliquee a la premiere 
masse. 
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5.2.1 Equations different ielles 

Les equations differentielles du mouvement s’ecrivent : 

mx i + (k + K ) xi + axi — Kx 2 = F 
—Kx i + mx 2 + (k + K ) :r 2 + ax 2 = 0 

5.2.2 Etude du regime permanent sinusoidal 

Solution permanente 

La solution generale de systeme d’equations differentielles est egale a la solution de la solu- 
tion du systeme homogene et d’une solution particuliere. La solution de l’equation homogene, 
en raison de l’amortissement, tend vers zero lorsque le temps augmente. Lorsque le regime 
permanent s’etablit, la solution devient egale a la solution permanente et s’ecrit alors : 

X\ = X\ cos (fit + (j)^ 
x 2 = X 2 cos (Qt + (j) 2 ) 

Pour calculer les amplitudes Ah et X 2 , ainsi que les phases d, et </> 2 , utilisons la methodes 
des nombres complexes. On peut, ainsi ecrire : 

x i = Re (X, e jnt ) a; 2 = Re (X 2 e jnt ) F = Re (F e jQt ) 

Dans ces expressions les amplitudes complexes sont definies par 

Xi = Xi e j ^ X 2 = X 2 e j ^ F = F 0 e j0 

Dans ce cas les equations differentielles se transforment en equations algebriques : 

f (k + K - mQ 2 + jaQ) X x - KX 2 = F 
\ -KX, +(k + K - mO 2 + jaQ) X 2 = 0 
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Amortissement negligeable 



Considerons d’abord le cas d’un amortissement suffisamment faible pour que Ton puisse 
considerer que a ~ 0. Le systeme d’equations different ielles s’ecrit alors 

J (k + K-mV 2 )X 1 ~KX 2 = F 

1 —KX.\ + (k + K - mV, 2 ) V 2 = 0 



Les solutions de ce systeme sont : 

*1 = 



2C 2 = 



F (9 2 a - V 2 ) 
m (c of — V 2 ) (c j 2 — O 2 ) 
KF 1 

m 2 (ujj - V 2 ) (a; 2 - V 2 ) 



x i • / k k + 2K . . 

Les pulsations uj\ = \ — et cu 2 = \ sont les pulsations propres calculees au cnapitre 

V m V m 

precedent. La valeur de la pulsation Q,\ est : 



v A = J k -^- 

V TO 

Les amplitudes des deplacements X 1 et X 2 sont alors donnees par 

I? 102 o2| 





F 


Vi = 


TO C 0 

KF 


v 2 = 


TO 2 



Les variations des amplitudes X\ et X 2 sont representees sur les figures ci-dessous 




Variation de Xi en fonct-ion de V Variation de X 2 en fonction de Q 



On remarque que le phenomene de resonance se produit pour Xi comme pour X 2 lorsque la 
pulsation d’excitation 0 est egale a l’une des pulsations propres oq ou ui 2 du systeme. L’amortis- 
sement etant tres faible, les amplitudes a la resonance sont tres importantes. Lorsque la pulsa- 
tion Q devient tres grande, ces amplitudes tendent vers zero. Enfin lorsque Q = ff i, l’amplitude 
Xi est egale a zero ; pour cette raison, la pulsation V A est appelee pulsation d’antiresonance. 
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5.3 Impedance 

Considerons le systeme a deux degres de liberte etudie dans le paragraphe precedent dans 
lequel nous supposons que l’amortissement est nul (a ~ 0). En regime stationnaire, on obtient 
pour l’amplitude complexe de la vitesse A , : 

X = : tt F 

— 1 J m(W -u\)~ 

On en deduit l’impedance d’entree : 

v F .m (O 2 — cuf) (O 2 — cu|) 

- E = 

Les figures ci-dessous donnent les variation de X et Z E en fonction de 0. On note le phe- 
nomene de resonance lorsque la pulsation d’excitation Q est egale a l’une des deux pulsations 
propres oj\ on A ces pulsations, le module de l’impedance d’entree est nul. Enfin, lorsque 
0 est egale a la pulsation d’antiresonance O i, la vitesse de la premiere masse est nulle et le 
module de l’impedance d’entree est infini. Lorsque O — > 00 , Z E ~ mil. 




5.4 Application 

Le phenomene d’antiresonance peut etre avantageusement utilise pour supprimer une vibra- 
tion resultant d’une resonance dans un systeme mecanique. 

Considerons le systeme a deux degres de liberte de la figure ci-dessous. 
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A 

Deux degres de liberte Variation de Xi en fonction de il. 

Les equations differentielles du mouvement s’ecrivent 

J rax i + axi + (k + K ) X\ — Kx 2 = F 
\ —Kx i + mx 2 + Kx 2 = 0 



n 



En regime permanent sinusoidal, on obtient 

F 0 



X i = 



*2 = 



D 2 - — 

“ M 



"> [° 4 + « 2 (br + b) - S] + V - § ) 

KF„ 1 



Mm [f)4 + £*> (M + f) - S] + j£fi (£*> - f ) 



Lorsque la pulsation de la force excitatrice est egale a cc.i = y/||, la masse m est immobile 

(Xi = 0). 

Si on choisit K et M telles que F = K (c’est-a-dire telles que cu 0 = &a), la masse m est 

immobile lorsque la pulsation excitatrice Q est egale a cu 0 = Dans ces conditions, 

l’ajout de M et K permet d’annuler la vibration de m a cette pulsation. Un tel dispositif 
constitue un ’’etouffeur” dynamique de vibrations. 
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Chapitre 6 

Generalites sur les phenomenes de 
propagation 



6.1 Propagation a une dimension 



6.1.1 Equation de propagation 

Dans les phenomenes vibratoires traites dans les chapitres precedents, nous nous sommes 
interesses a des irhenomenes ou des grandeurs physiques qui dependaient d’une seule variable, 
le temps. Nous allons maintenant examiner t-oute une une serie de phenomenes qui sont decrits 
par une fonction qui depend a la fois du temps t et d’une variable d’espace , x par exemple. 

Ces phenomenes sont regis par une equation aux derivees partielles, appelee equation d’onde 
ou equation de propagation a une dimension de la forme : 



d 2 s 1 d 2 s 

dx 2 V 2 dt 2 



( 6 . 1 ) 



dans laquelle V est une grandeur physique qui a les dimensions d’une vitesse et sera appelee 
dans la suite vitesse de propagation. 



6.1.2 Solution de l’equation de propagation 

Methode de D’Alembert 

Pour resoudre Pequation des ondes a une dimension, operons le changement de variable 
suivant : 



V = (6.2) 

e = (6.3) 

Calculons les derivees partielles par rapport a t et x, en fonction des derivees partielles par 
rapport a 77 et £. 
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Sachant que : 



on obtient 



dr/ <9£ dr i 1 

dt^di = et que 0 x = ~£h: = V 



ds ds dr] ds d £ ds ds 

dt dr] dt <9£ dt dr] <9£ 

ds ds dr/ ds <9£ 1 ds ds 

dx dr] dx d£ dx V dr] <9£ 

En tenant compte de ces resultats et sachant que 

d 2 s d 2 s 



dr/d £ d£dr] 



on obtient : 



rp 1 „ pjl ~ f\ 2 „ Cj 2 ~ 

W = ~^~ 2 d^ + d^ ^ 

d 2 s _ 1 \d 2 s o d 2 s <9 2 sl 

dx 2 ~ C 2 \g~2 ~ 2 d^di + d^\ ^ ' ’ 

En remplagant dans l’equation d’onde |)Sf et |^f par les expressions ci-dessus, on obtient 
l’equation d’onde exprimee en fonction des derivees partielles par rapport aux variables r] et £ : 

d 2 s 



dr/d £ 

Cette derniere equation peut s’ecrire 

d ds 
dr] <9£ 

Un integration par rapport a 7/ donne : 



ds 

rr fM 



( 6 . 10 ) 



( 6 . 11 ) 



( 6 . 12 ) 



ou / ( 7 /) est une fonction qui ne depend que de // (et pas de £). Enfin une integration par rapport 
a £ donne : 

s(,,e) = F(i|)+e(5) (6.13) 

ou F (//) , qui ne depend que de 7 /, est une primitive de / (77). La fonction G (£) est une fonction 
qui ne depend que de £. En revenant aux variables x et t. on obtient la solution generate de 
l’equation des ondes a une dimension : 



s ( x , t ) = F (t - + G (t f ^ 



(6.14) 



Les fonctions F (t — |f) et G (t + |t) sont des fonctions dont la nature est fixee par les 
conditions aux frontieres imposees a la solution s (x,t). 
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Proprietes des solutions particulieres F (t — y^j et G (t + y) 

Proprietes de F (t — y) On etudie le cas de la solution particuliere F (t — y). Pour cela 
on suppose que les conditions aux frontieres sont telles que G (t + y) est constamment nulle. 
On considere a l’instant t\ un point d’abscisse x\ . La valeur de la fonction s en ce point et a cet 
instant est s (xi,ti). On recherche a un instant t 2 posterieur a t\ (t 2 > t\) la position x 2 d’un 
point pour lequel la valeur de s est la meme que la valeur qu’elle avait en x\ a l’instant t \ . Ce 
probleme est formule par l’egalite suivante : 





S (Xyti) 


= s(x 2 ,t 2 ) 


(6.15) 


Ce qui se traduit par 

F 

Cette equation est satisfaite si 






(6.16) 








(6.17) 


D’ou la valeur de x 2 : 


x 2 = xi + V (t 2 - ti) 


(6.18) 


Comme t 2 > t\, x 2 est superieure a X\ et ces deux points sont dist-ants de 






X 2 —Xi = 


~-V(t 2 - ti) 


(6.19) 



F (t — y) correspond a une onde se propageant- dans le sens des x croissants (Voir la figure 
ci-dessous). F (t — y) est appelee onde progressive et cette expression constituera dans la suite 
la definition d’une onde progressive. 



t=t, 




Direction de 
propagation 






X2 




< *• 

x 2 -xi=V(t 2 -ti) 
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Proprietes de G (f + f-) On etudie le cas de la solution particuliere G (t + y). Pour cela 
on suppose que les conditions aux frontieres sont telles F (t — est constamment nulle. On 
considere a l’instant t \ un point d’abscisse x \ . La valeur de la fonction s en ce point et a cet 
instant est s (x On recherche a un instant t 2 posterieur a t\ (f 2 > f i ) la position a : 2 d’un 
point pour lequel la valeur de s est la meme que la valeur en x± a l’instant t\ . Ce probleme est 
formule par l’egalite suivante : 





s(xi,t 1 ) = s(x 2 ,t 2 ) 


(6.20) 


Ce qui se traduit par 


( fl + v) = G (* 2 + v) 




O 


(6.21) 


Cette equation est satisfaite si 


X\ x 2 

tl + V = t2 + V 


(6.22) 


D’ou la valeur de x 2 : 


x 2 = xi-V (t 2 - fi) 


(6.23) 


Comme f 2 > t\, x 2 est inferieure a X\. Ces deux points sont distants de 






x 1 -x 2 = V (t 2 - fi) 


(6.24) 



G (t + y) correspond a une onde se propageant dans le sens des x decroissants (Voir la 
figure ci-dessous). G (t + correspond a une progressive se propageant- dans le sens des x 
decroissant. 





xi-x 2 =V(t 2 -ti) 



6.1.3 Onde progressive sinusoidale 

On considere une onde progressive se propageant dans la direction de l’axe des x. telle que 
le point d’abscisse x = 0 est soumis a une vibration sinusoidale de la forme 

s (x = 0, t) = So cos (cut) (6.25) 
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Le point se trouvant a l’abscisse x > 0 aura la meme vibration que celle du point x = 0 
mais avec un retard egal a y : 



s (. x , t) = So cos 



u 



X 

'~v 



(6.26) 



Cette expression constit.ue la definition d’une onde progressive sinusoi'dale (ou harmonique) ; 
elle peut et-re ecrite sous la forme : 



s (. x , t) = So cos [cut — (f) (x) 



(6.27) 



ou 4> (x) = yx represente le dephasage lie au temps de propagation y. On dit que <i> (x) 
represente le dephasage du a la propagation. L’onde progressive sinusoi'dale s’ecrit sous la forme 
suivante qui permet de mettre en evidence la double periodicite (dans le temps et dans l’espace) : 



s (x, t ) = S 0 cos 





(6.28) 



La quantite T = — est la periode temporelle t-andis que la quantite A 
d’onde qui constitue la periode spatiale. On peut verifier aisement que : 



VT est la longueur 



s(x,t + nT ) = s(x,t ) 
s (x + nX) = s(x,t ) 



ou n est un nombre entier. 

L’onde progressive s’ecrit souvent : 

s (x, t) = So cos [ cut — kx] 

ou k = y = est- appele le module du vect-eur d’onde qui s’exprime en rn 1 . 

On utilise tres souvent la notation complexe d’une onde progressive sinusoi'dale : 

s (x, t ) = So e i{ut ~ kx) 
s (x, t) = Se iujt 



(6.29) 

(6.30) 



(6.31) 



(6.32) 

(6.33) 



ou S_ = S (J e~ lkx represent-e l’amplit-ude complexe de l’onde progressive sinusoi'dale. Le 
module So de S_ est l’amplitude de l’onde tandis que son argument —kx represent-e le dephasage 
du a la propagation. 



6.1.4 Superposition de deux ondes progressives sinusoi'dales 

Cas de deux ondes de meme frequence se propageant dans le meme sens 

Considerons deux ondes de meme frequence et- de meme direction de propagation, d’ampli- 
t-udes respectives ,S'i et ,S' 2 , et de phases respectives d>\ et- </> 2 . L’onde resultante sera alors : 

s {x, t) = Si + s 2 e j{^t-kx+4> 1 ) = s e j(ut-kx+4 > ) (6.34) 
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ou encore en notation reelle : 



s (x, t) = S cos (cut — kx + 4 i>) 



(6.35) 



avec 



et 



S = 




+ <Sf + 2 S 1 S 2 COS (</>! 




(6.36) 



0 = Arctg 



f sin (0i) + S 2 sin (0 2 ) 
\Si cos (0 : ) + 5*2 cos (0 2 ) 



(6.37) 



La superposition de deux ondes harmoniques de meme frequence, et qui se propagent dans la 
meme direction, donne une autre onde harmonique progressive de meme frequence, d’amplitude 
S et de phase 0. 



Cas de deux ondes de meme frequence se propageant dans des sens opposes 

Si par contre, on superpose deux ondes harmoniques de meme frequence mais se propageant- 
dans des sens opposes , le resultat est tout autre. En effet, dans ce cas : 



s (x, t ) = e + S 2 eP^+kx+h) = ^s^e~ ikx + S 2 e j ^e +jkx ] e jut (6.38) 

et on ne plus ecrire l’onde resultante sous la forme d’une onde progressive simple. Un cas 
particulier important se produit quand les deux amplitudes sont identiques. Si on note : 

S 1 = S 2 = So (6.39) 

on a : 

5 ( x , t) = 2 S 0 cos (kx + ^ (6.40) 

et done en notation reelle : 

02 



Ce mode de vibration est tres different d’une onde progressive puisque tous les points x de 
la corde vibrent en phase avec des amplitudes differentes. En particulier, il existe une serie de 
points : 




s (x, t ) = 2Sq cos ( kx + 



cos cut 



+ 



avec 





(6.42) 



http:/ / dj elouah.ifrance.com 






6.1 Propagation a une dimension 



51 



n = 0, ±1, ±2, 

ou 1’ amplitude de vibration est constamment nulle. On dit dans ce cas que l’onde est sta- 
tionnaire et que les points x n sont les noeuds de l’onde. Entre chaque paire de noeuds existe un 
ventre ou l’amplitude de vibration est maximum et egale a 2 So. On note aussi que l’intervalle 
entre deux noeuds est egal a une demi-longueur d’onde A/2. 

6.1.5 Vitesse de phase 

Considerons une onde progressive sinusoidale qui se propage dans le sens des x croissant. 
Un point d’abscisse x possede, a l’instant. f, l’elongat.ion : 

s ( x , t ) = S 0 cos ( cut — kx ) (6.43) 

Entre l’instant t et t + At l’onde progresse d’une quantite Ax. A l’instant. t + At, le point 
d’abscisse x + Ax possede la meme elongation que celle que possedait- le point d’abscisse x a 
l’instant. anterieur t. Ceci se traduit par l’egalite : 

s(x,t) = s(x + Ax,t + At) (6.44) 

S 0 cos (cut — kx) = S 0 cos [cu (t + At) — k(x + Acc)] (6.45) 

Cette egalite est satisfaite si les phases instant-anees sont egales : 

cut — kx = cu (t + At) — k (x + Ax) (6.46) 

Soit encore 

cu At = k Ax (6.47) 

On definit- la vitesse de phase V/ = df qui s’exprime en fonction de cu et k apr : 

U = j (6.48) 

Si la vitesse de phase ne depend pas de cu, le milieu est dit non dispersif. Dans le cas contraire 
il est dit dispersif. 

La figure ci-dessous permet d’illustrer la notion de vitesse de phase en considerant deux 
representations a des instants diflerents d’une corde parcourue par une onde . La courbe continue 
represente l’ensemble des points de la corde a l’instant t. Le point de la corde d’abscisse x est 
represente par le point blanc, t-andis que le point d’abscisse x + Ax est represente par le point 
noir. On constate qu’ent-re les instants t et t + At chacun de ces point suit une trajectoire 
rectiligne et le deplacement du point noir a l’instant t I At est egal au deplacement du point 
blanc a l’instant t. En particulier la crete de la corde, correspondant a une valeur particuliere 
de la phase inst-antanee, semble se deplacer dans le sens de propagation de l’onde avec la vitesse 
de V/, mais la trajectoire de chaque point materiel est une trajectoire rectiligne perpendiculaire 
a la direction de propagation. 
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6.1.6 Vitesse de groupe 

La vitesse de phase n’est pas necessairement la vitesse que Ton observe lorsqu’on analyse 
un mouvement ondulat-oire. En general une onde n’est pas parfaitement sinusoi'dale mais a une 
duree limitee et se presente sous la forme d’un train d’onde appele communement ’’pulse” ou 
’’groupe” qui se propage avec une vitesse Vq appelee vitesse de groupe. Cette onde sous la 
forme d’un pulse cont.ient plusieurs frequences. Si la vitesse de phase est independante de la 
frequence (Milieu non dispersif) alors toutes les frequences qui constituent le pulse se propagent 
a la meme vitesse et le pulse se propage avec une vitesse de groupe egale a la vitesse de phase. 
Mais si le milieu est dispersif (i.e la vitesse de phase depend de la frequence), alors le pulse se 
propage avec une vitesse de groupe differente de la vitesse de phase. 

Pour illustrer ce phenomene, considerons une onde constitute de deux ondes de frequence 
differente et de meme amplitude. En x = 0, cette onde s’ecrit par exemple sous la forme : 



S (0, t) = So COS (iUlt) + So COS (o>2 1) 


(6.49) 


Cette onde peut s’ecrire encore : 




s (0, t) = 2S 0 cos (ust) cos (cut) 


(6.50) 


ou 

UJ2-UJ! ^ 0J 2 +CUi 

uj b = 2 et u = 2 


(6.51) 



Si oji est voisine de gv 2 , la vibration resultante se presente sous la forme d’une sinusoi'de 
de pulsation oj dont l’amplitude est modulee par un battement de pulsation lu b (Modulation 
d’amplitude). 
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En un point x > 0. l’onde obtenue resulte de la superposition de ces deux ondes qui se sont 
propagees a des vitesses differentes car le milieu de propagation est suppose dispersif : 

s (x, t ) = S 0 cos (uit. — k\x) + S 0 cos (oj 2 t — k 2 x) (6.52) 

s (x, t ) pent s’ecrire : 



s ( x , t) = 2^0 cos (c ust — ksx) cos (cut — kx ) 
Dans cette expression : 

k 2 — k\ k 2 + ki 

k B = — - — et k = — - — 



(6.53) 



(6.54) 



L’amplitude du battement se propage a une vitesse qui est la vitesse de groupe definie par 
la relation : 

UJ B CU2-CU1 du . . 

G “ k B - k 2 -h ~ dk ^ ^ 

Comme cu 2 est peu diflerente de oq . la vitesse de groupe est definie par : 



V G = - 
G dk 



(6.56) 



Tandis que la sinusoi'de contenue a l’interieur du battement se propage a la vitesse de phase : 



(6.57) 




Les fleches verticales noires correspondent au maximum des battements qui se propagcnt a la 
vitesse de groupe. Les fleches verticales blanches correspondent au maximum des vibrations 

qui se propagent a la vitesse de phase. 



http:/ / dj elouah.ifrance.com 







6.2 Propagation en trois dimensions 



54 



6.1.7 Onde Vectorielle 



Dans ce qui precede, la quantite s ( x , t) represente une grandeur scalaire, mais certains 
phenomenes decrits par des vecteurs conduisent a des equations similaires : 

d 2 A 1 d 2 A 



dx 2 V 2 dt 2 



Le vecteur A, defini dans un milieu trois dimensions, a trois composantes A x . A y . A z et l’ex- 
pression ci-dessus signifie que chacune de ces composantes satisfait individuellement l’equation 
de propagation : 

d 2 A x 1_ d 2 A x _ n 

dx 2 V 2 dt 2 

= 0 ( 6 ' 59 ) 

&*A Z j_ 9*A Z _ n 

dx 2 v 2 at 2 u 

Chacune de ces composantes A x , A y , A z se propage en t — ^ et t + 



(6.58) 



6.2 Propagation en trois dimensions 



6.2.1 Equation de propagation 

Dans un systeme de coordonnees cartesiennes, l’equation de propagation en trois dimensions 
s’ecrit sous la forme : 

d 2 s d 2 s d 2 s 1 d 2 s . . 

dx* + £hf + d^ 2 ~ V 2 dt 2 = ' ' ^ 

On definit le laplacien scalaire de s par l’expression ci-dessous : 



_ d 2 s d 2 s d 2 s 

dx 2 dy 2 dx 2 

et l’equation des ondes s’ecrit sous la forme condensee : 

A 1 d 2 s n 

As ~v 2 di 2 ~° 



(6.61) 



(6.62) 



s est fonction du temps mais egalement des cordonnees du point M ou la fonction s doit 
etre calculee. Si Ton appelle le vecteur position f = OM . la quantite s depend du temps et du 
vecteur position f ; on ecrit s ( r,t ). 



6.2.2 Onde plane progressive sinusoidale 

Definition 

L’onde progressive sinusoidale (ou harmonique), se propageant dans une direction donnee 
par un vecteur unitaire u est definie par : 

s (r, t) = Sq cos (cut (6.63) 
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ou le vecteur k = k u est appele le vecteur d’onde. Si les composantes du vecteur k sont k x . 
k y et k z . alors l’onde plane est definie par 



s ( x , y, z , t ) = S 0 cos ( cut — k x x — k y y — k z z) (6.64) 

On peut utiliser la notation complexe pour representer l’onde plane progressive sinu- 
soi'dale qui s’ecrit dans ce cas sous la forme : 

s(f,t) = S 0 (6.65) 



Relation de dispersion 

En remplagant s(r,t ) par son expression dans l’equation de propagation, on obtient la 
relation k = k (cu) pour que l’onde plane definie ci-dessus constit.ue une solution particuliere de 
l’equation d’onde. Cette relation est appelee la relation de dispersion et elle s’ecrit : 

k = — (6.66) 



Surface d’onde 

On appelle surface d’onde ou surface equiphase, l’ensemble des points de l’espace pour 
lesquels, au meme instant, s (r, t) a la meme valeur. Recherchons la surface d’onde passant par 
un point M 0 a un instant t ; cette surface est l’ensemble des points M de l’espace pour lesquels 
l’egalite suivante est sat.isfaite : 



s(r,t) = s(f 0 ,t) 


(6.67) 


Cette egalite se traduit par : 




S 0 = S 0 


(6.68) 


Cette egalite est satisfaite si 

k ■ f = k ■ f 0 


(6.69) 


Tenant compte des proprietes du produit scalaire, on obtient 




k ■ (f — fo) = 0 


(6.70) 



La surface d’onde passant, a l’instant t, par le point M 0 est l’ensemble des points M satis- 
faisant. l’equation ci-dessus. Cette surface est un plan passant par M 0 , et perpendiculaire a la 
direction du vecteur d’onde, done a la direction de propagation. L’onde est dite plane. 
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II existe d’autres types d’ondes definis par les surfaces d’onde respectives : par exemple les 
ondes spheriques pour lesquelles les surfaces d’onde sont des spheres ou les ondes cylindriques 
pour lesquelles les surfaces d’onde sont des cylindres. 

Polarisation 

Dans le cas d’une onde plane progressive sinusoi'dale representee par une quantite vectorielle 
A ( r,t ), cette quantite peut avoir differentes orientations par rapport aux surfaces d’ondes : 

1. A est constamment perpendiculaire a la surface d’onde, ou de maniere equivalente paral- 
lel a la direction de propagation : l’onde est dite longitudinale. 

2. A est contenu dans la surface d’onde, ou de maniere equivalente perpendiculaire a la di- 
rection de propagation : l’onde est dite transversale. Dans ce cas, Pextremite du champ 
vectoriel A peut decrire une trajectoire rectiligne : l’onde transversale est dite a polarisa- 
tion rectiligne. Elle peut decrire une trajectoire circulaire (onde transversale a polarisation 
circulaire), ou une trajectoire elliptique (onde transversale a polarisation elliptique). 
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Chapitre 7 

Cordes vibrantes 



7.1 Equation des ondes 

Considerons une corde tendue, rectiligne selon la coordonnee x, et de longueur infinie. Nous 
allons etudier la propagation d’un faible ebranlement le long de la corde. Supposons que cet 
ebranlement se produise suivant l’axe Oy. 

Etudions b equation du mouvement de cette corde. Nous denoterons par T la tension a 
laquelle est soumise la corde. On considere en un point d’abscisse x un segment tres court de 
cet-te corde, de longueur Ax. La masse Am du segment est donnee par : 

Am = jiAx 

ou n est la densite lineique de masse de la corde, c’est-a-dire la masse par unite de longueur 
qui s’exprime en kg/m. 




Corde vibrant transversalement 



Dans une situation hors equilibre, le segment n’est plus droit, il presente une courbure. Nous 
considerons des mouvements d’oscillation de la corde de petite amplitude 
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u (. x , t) = u (. x , t) e y 

si bien que nous pouvons faire 1’ approximation : 



sin (9) \ x = tan ( 9 ) 



du 

dx 



X 



sin (0)L+A* = tan ( 0 )L+A* = 



du 

dx 



x-\-Ax 



Cette approximation neglige aussi l’allongement du segment, et considere done la tension T 
comme constants. La force appliquee sur le segment dans la direction y est la resultante de la 
force appliquee au point x, qui est une force appliquee vers le bas et egale en module a 



F (. x , t) = T sin 



= T tan 



du 
x = dx 



et de la force appliquee au point x+Ax qui est vers le haut et egale a 



F (x + Ax, t ) 



T sinWL+A* - T tan ( e )\ x+ Ax 



du 

dx 



x-\-Ax 



La force t.otale dans la direction y est done : 



R 



F (x + Ax, t) — F ( x , t) =T 



du 


du 




dx 


x+Ax dx 


X- 



d^u 

dt 2 



Ax 



Nous pouvons appliquer maintenant la loi fondamentale de la dynamique au segment Ax. La 
force dans la direction y doit etre egale au produit de la masse Am du segment par l’acceleration 
de celui-ci. Done : 



R = Am 
d 2 u 






dt 2 



= T 



d 2 u 

~dF 

d 2 u 



dx 2 



Si on definit V = J ^ qui a la dimension d’une vitesse, on constate que : 



d 2 u 1 d 2 u 



= 0 



dx 2 V 2 dt 2 

qui est l’equation d’ondes de la corde. V qui est la vitesse de propagation de cette onde. 
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7.2 Ondes progressives harmoniques 

7.2.1 Definition 

Une onde progressive harmonique se propageant selon Ox est definie par : 

u (x, t ) = Uq cos (cut — kx ) 
ou encore en notation complexe 



u (x, t ) = U 0 e j(ut ~ kx) 

ou /; = | = j est le module du vecteur d’onde, A et-ant la longueur d ’onde. 

7.2.2 Force en un point 

On appelle force en un point, la projection selon Oy de la force exercee, en ce point, par la 
partie gauche de la corde sur la partie droite : 

f = - t 7T 

OX 

Dans le cas d’une onde progressive sinusoi'dale, cet-te relation devient : 

F(x,t) = jkT U 0 e j{ujt ~ kx) 

La vitesse de particules s’ecrit : 

I 

u{x,t) = — =juU Q 

On constate que pour une onde progressive la vitesse de particules u est en phase avec la 
force F. 



7.2.3 Impedance 

On appelle impedance en un point le rapport de l’amplitude complexe de la force a l’am- 
pliude complexe de la vitesse de particule 

Z (' X ) = Y 

Dans le cas d’une onde progressive, on obtient : 

Z (. x ) = i-iV = \J i-iT 

La quantite \f fi I definit 1’ impedance caracteristique de la corde 

Z c = \fyT = fj,V 
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On obtient une propriety de l’onde progressive plane 

Z (. x ) = Z c Wx 

7.3 Oscillations libres d’une corde de longueur finie 




Corde de longueur L fixee aux extremites 



Considerons une corde de longueur L fixe aux points x = 0 cl x = L. Recherchons une 
solution de l’equation d’onde sous la forme : 



u{x,t) =g(x)f(t ) 

En remplagant dans l’equation de propagation, on obtient : 

1 d 2 g 1 1 cl 2 f 

= V*f Id? 

Le membre de gauche de cette equation ne depend que de x, t-andis que le membre de gauche 
ne depend que de t. Ces deux expressions sont done egales a une constante qui doit et-re un 
nombre reel negatif que nous posons egal a —k 2 car la solution ne doit pas tendre vers l’infini 
lorsque t tend vers l’infini. Posons uj = kV . On en deduit que : 



d 2 g __ 

£ = 

dt - 2 



-k 2 g 

~0J 2 f 



Les solutions de ces deux equations differentielles sont de la forme : 



/ = A cos (cut) + B sin (cut) 
g = C cos (kx) + D sin (kx) 
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La solution de l’equation d’onde peut alors s’ecrire sous la forme : 

u ( x , t ) = [A cos (cut) + B sin (cut)] [C cos (kx) + D sin (kx)] 
Tenant compte des conditions aux limites 

u (0, t) = 0 
u (L, t) = 0 

on obtient 



C = 0 

k = rij^ o n = 0, 1, 2, 

La solution de l’equation d’onde qui satisfait ces conditions aux limites est done une somme 
d’une infinites de termes : 



avec 



OO 

u ( x , t) = [a n cos (cu n t) + b n sin (cu n t)] sin (k n x) 

n = 0 



k n = n~r et cu n = k n V = n— - 

1j 1j 

Les cu n sont les pulsations propres. Les coefficients a n et- b n sont determines par les conditions 
initiates du mouvement. Supposons qu’a t = 0 nous imposions a la corde une cert-aine forme 
initiate u(x, 0) = u 0 (x ) et une vitesse initiale 



u (0, t) = v 0 (t) 

Dans ce cas nous aurons les conditions initiales suivantes : 



OO 

U 0 (x) = a n sin (Kx) 
n = 0 
oo 

^0 (x) = -u n b n sin (k n x) 

n = 0 

On doit inverser ces equations pour obtenir les coefficients a n et b n . La methode de Fourier 
consiste a les multiplier par sin ( k m x) et les integrer ent-re 0 et L. Si on utilise les integrales : 



on obtient 




dx 



0 si m 7^ n 

L 

— si m = n 

Zj 



Cl n 



b 



n 



— / u 0 [x) sm 

L J 0 

o rn 



cu„L 



n 0 (x 




sin 




dx 
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7.4 Reflexion et transmission 

7.4.1 Reflexion et transmission entre deux cordes semi-infinies 




incidente reflechie t 



transmise 




Z 1 = ^[m{T z 2 = r 



Reflexion transmission dans deux cordes semi-infinies 



Soit deux cordes de longueur semi-infinie, reliees en x = 0. Leurs masses lineiques sont 
respectivement // 1 et // 2 . Lorsqu’une onde venant de — oo se propage vers x = 0 dans la premiere 
corde, elle donne naissance au point de jonction, x = 0. a une onde reflechie et une onde 
transmise. L’ecriture de la continuite du deplacement et de la force en x = 0 permet d’obtenir 
le coefficient de reflexion 11, , et le coefficient de transmission T u definis respectivement par : 



r 

” ” 

T — — 

u ~ Ui 

ou Ui , Ur et Ur sont les amplitudes des deplacements associes respectivement a l’onde 
incidente, l’onde reflechie et l’onde transmise. On en deduit : 



R 



U 



Z\ — Zi2 
Z\ + Z2 



T 



2 Z 1 

Z\ + Z 2 



7.4.2 Reflexion sur une impedance quelconque 
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z c - 



Corde semi-infinie t-erminee par une impedance Z^ 



Soit- une corde de longueur semi-infinie, de masse lineique //. tendue horizontalement avec une 
tension T et t-erminee en x = 0 par une impedance mecanique Zj. Lorsqu’une onde harmonique 
se propage dans la corde de — oo vers x = 0. elle subit une reflexion en ce point. Sachant que le 
deplacement de particules s’ecrit : 



u (x, t ) = Ui e j ^~ kx) + U R e j ^ t+kx) 

on en deduit la vitesse de particules et la force en un point d’abscisse x 



u = ^ = joj [Ui e^ ut k U + Ur e^ ut + kx )~^ 

F = -Tg = -jkT [Ui UUi-kx) _ Ur e j{wt+kx )] 



En x = 0, les conditions aux limit-es s’ecrivent 



Zt — 



mt) 

u{ 0, t ) 



On en deduit le coefficient de reflexion R u en fonct-ion de l’impedance caracteristique Z c et 
de l’impedance Zr placee a l’extremite de la corde : 



R 



U 



Ur 

Ui 



Z c — Z T 
Zc + Zr 
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Chapitre 8 

Ondes acoustiques dans les fluides 

8.1 Introduction 



Les ondes acoustiques sont des ondes elastiques qui se propagent dans les fluides (gaz ou 
liquides). II est done possible d’obtenir l’equation d’onde qui regit la propagation des ondes 
planes dans un fluide par la meme demarche que celle que nous avons utilisee pour etablir 
l’equation de propagation des ondes transversales dans une corde. 

Dans la suite, nous utiliserons les symboles suivants pour etudier l’onde acoustique qui se 
propage suivant l’axe des x : 

x : coordonnee a l’equilibre d’une particule du milieu. 

u x : composante suivant l’axe des x du deplacement de particule par rapport a la position 
d’equilibre. 

p 0 : masse volumique du fluide a l’equilibre 
P : pression instantanee en un point quelconque 
P 0 : pression a l’equilibre 

p = P — Pq : surpression ou pression acoustique 
c : vitesse de propagation de l’onde 

On entend par particule, un element de volume contenant des millions de molecules de 
telle sorte qu’il puisse etre considere comme continu, mais toutefois suffisamment petit pour 
que les grandeurs acoustiques comme la pression, la masse volumique et la vitesse de particule 
puissent etre considerees comme const.antes dans cet element de volume. Dans ce qui suit, nous 
negligerons les effets de la gravitation de telle sorte que Iq et p 0 sont uniformes dans tout le 
milieu. On suppose d’autre part que le milieu est homogene, isotrope et parfaitement elastique, 
e’est-a-dire non dissipat.if. 



8.2 Equation d’onde 

Considerons le cas d’une onde plane emise dans un fluide par une membrane vibrante plane. 
Lorsque celle-ci est au repos, la pression dans le fluide est uniforme et egale a Pq. En se deplagant, 
par exemple dans le sens des x positifs, la membrane comprime la couche de fluide adjacente. 



http:/ / dj elouah.ifrance.com 



64 





8.2 Equation d’onde 



65 



Cette situation est instable : le fluide se defend en comprimant a son tout la tranche voisine. 
L’onde progresse ainsi de proche en proche par une succession de compressions et de detentes. 



P(X) . 


1 

1 

1 

1 

1 

1 

1 


1 

.P(x+dx) ! 


n 

s 

] 




i 

1 

i 


X 




x+dx 






► 


► 






u x 0) 


u x (x+dx) 





◄ ► 

dx 

Propagation acoustique d’une onde acoustique 



Soit une tranche de fluide de petite epaisseur Ax sit-uee a l’abscisse x. Lorsque la pertur- 
bation att-eint ce point, les forces agissant sur cette tranche ne s’equilibrent plus et elle se met 
en mouvement. Soit u x ( x , t) le deplacement a l’instant t du plan d’abscisse x. Soit F x (x, t ) et 
F x (x + Ax, t ) les forces agissant sur la tranche de fluide respectivement en x et x - 1 Ax. Ces 
forces s’expriment par : 



F x (x, t) = SP (x, t ) 

F x {x + Ax, t) = — SP {x + Ax, t) 

La result-ante de ces deux forces est : 

A F x = F x (x, t ) + F x (x + Ax, t ) 

A F x = —S [P (x + Ax, t) — P (x, f)] 

En faisant- un developpement en serie de Taylor au premier ordre de P (x, t). on obtient : 

dP 

P(x + Ax, t) = P (x, t) + — — Ax 

ox 

D’ou : 



dP 

A F x = -S — 

Ox 

Comme P = Pq + p, la force resultant s’exprime par : 



AF x = ~S^Ax 
Ox 
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Sous Taction de cet-te force, la tranche de fluide subit une acceleration et en ecrivant la 
relation fondamentale de la dynamique, on obtient : 



Am 



d 2 u x 

~W 



„dp 

dx 



Ax 



Le fluide etant compressible, le deplacement du plan d’abscisse x + Ax est different du 
deplacement du plan d’abscisse x et il vaut u x (x + Ax. t). De nouveau un developpement en 
serie de Taylor au premier ordre permet d’ecrire : 

, . . . . du x 

u x [x + Ax, t) = u x (x, t) + -q— Ax 

Pour prendre en compte la compressibilite du fluide, calculons la dilatation volumique subies 
par la tranche de fluide . Soit Ary = S Ax, le volume a Tequilibre et soit Av, le volume en 
cours de mouvement, avec : 



Au 

Au 

Au 

Au 



S [x + Ax + u x (x + Ax) — x — u x (a?)] 
S [Arc + u x (x + Arc) — u x (rc)] 



S 
Au 0 



Arc + Ax 



dui 

dx 



dx 

Aun 



On en deduit la dilatation volumique 



Au — Ary 
Ary 



du 

Rappelons que pour un fluide compressible, la surpression p est reliee a la dilatation volu- 
mique 0 par la relation 

p = — k 6 

ou k est le module de compressibilite. On obtient ainsi : 



p = — n 



du x 

dx 



Remarque : On utilise souvent. le coefficient de compressibilite x 
Les deux equations 



1 

K 
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d 2 u x dp 

Po ~dP = ~&x 

du x 

P = -k — 
ox 

constituent les deux equations fondamentales de l’acoustique. 

En remplagant dans la premiere equation p par son expression tiree de la seconde equation 
on obtient P equation de propagation : 



d 2 p 1 d 2 p 

dx 2 V 2 dt 2 

ou V represente la vitesse de propagation definie par 




V do VdoX 



8.3 Vitesse du son 

Le phenomene de propagation et-ant un processus adiabatique, la relation liant la pression 
et le volume est 

Pv Y = cons tan te 

En calculant- la differentielle, on obtient : 

v 7 dP + P v 1 ^ 1 dv = 0 



Si l’on considere que dP represente la varaiat-ion de pression au voisinage de la pression a 
Pequilibre i y. on obtient : 

v^ p + qP 0 Av = 0 

d’ou : 



p = -7P0 



Av 

Vo 



En tenant compte de la definition du module de compressibilite, on obtient : 

1 



« = - = 7^o 
X 



D’ou la vitesse du son dans un fluide : 



V = 



IjPo 

do 



Exemple : Dans Pair, dans les conditions normales T = 20° 6' et P 0 = 10 5 iV ■ m 2 , 7 = 1.4 
et p 0 = 1.29 kg ■ m~ 3 , on en deduit V ~ 330 m ■ s^ 1 
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La valeur de la pression a l’equilibre depend fortement de la temperature. Pour une mole 
de gaz par fait, on a : 



d’ou 



PqVq — RT 



Po 



RT 

Vo 



et 



V = 



IjRT 

Po v o 



Le produit p 0 vo represente la masse molaire M du gaz ; d’ou : 



1 / = 




Dans un gaz parfait, la vitesse de propagation du son est proportionnelle a la racine carree 
de la temperature mesuree en °K. 



8.4 Onde progressive sinusoi'dale 

8.4.1 Definition 

Une onde acoustique sinusoi'dale, s’ecrit : 

p ( x , t) = po cos 



, x \ 



On definit- le module du vecteur d’onde k par 



k~- 

k ~v 



d’ou 



p ( x , t) = po cos (c ot — kx) 

En notation complexe, l’onde progressive sinusoi'dale s’ecrit 

p ( x , t) = p 0 e j{ujt ~ kx) 

La relation liant la pression acoustique et la compressibilite, a savoir 

du 

P = ~ k 7T 
ax 
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permet d’ecrire 



u(x,t ) = 



1 

K 



p ( x , t ) dx 



u 



(x,t) = -- f poe^^-^dx 



u(x,t ) = -E2_ei("t-*a0 
jkn 

u(x,t) = 

jixp 0 V 

La derivation de cette derniere expression par rapport au temps permet d’obtenir la vitesse 
de particules : 



*(!,«) = ^ = T Lpoe« ( " , - to > 



dt p 0 V ‘ 

On constate que pour une onde progressive la vitesse de particules est en phase avec la 
pression acoustique. 



8.4.2 Impedance acoustique 

On appelle impedance acoustique en un point le rapport de Pamplitude complexe de la 
pression a l’ampliude complexe de la vitesse de particule 



Dans le cas d’une onde progressive, on obtient : 

Z(x) = p 0 V 

Le produit p 0 V definit- Pimpedance acoustique caracteristique du fluide 

= p 0 V 

On obtient une propriety de l’onde plane progressive : 

Z (. x ) = Z c Wx 

8.4.3 Energie acoustique 

Densite d’energie cinetique 

Soit un petit element de volume v 0 dont le depelacement est u (x, t) et dont la vitesse est 
u(x,t ) ; il possede une energie cinetique 

1 . o 

E c = t;Po v ° u 
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TO 



On definit l’energie cinetique par unite de volume ou densite d’energie cinetique 

E r 1 



£c — — o P 0 U 

Vo 2 



•2 



Densite d’energie potentielle 



Soit, un petit element de volume vq. Sous l’act-ion de la surpression p, cet element se comprime 
ou se dilate en raison de la compressibilite du fluide. L’energie potentielle emmagasinee est egale 
au travail fourni par la pression pour comprimer ou dilater le volume Vq : 



Sachant que 
on en deduit que : 

d’ou : 



E p = —pclv 



p = — K ■ 



v - vq 

Vo 



dv = Vo dp 

K 



Vo 



Ep = — pdp 

K Jo 

Ep = 

On en deduit la densite d’energie potentielle : 

g Ep 1 2 

L 'P n " 

Vq 2k 

Densite d’energie 



La densite d’energie est egale a la somme de la densite d’energie cinetique et de la densite 
d’energie potentielle 



£ — £ c + £ p 

C 1 • 2 I 1 2 

f = 2 P ° U + 2^ P 

Dans le cas particular d’une onde plane progressive sinusoi'dale, 

1 



£ = £ = 
Lsr i-"n 



2p„V’- 



po cos" (cut — kx) 
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et 



£ = 



1 



P,V 2 



p 0 cos 2 (cut — kx) 



On definit la valeur moyenne temporelle de la densite d’energie comme 

r t+T 



1 f 

{£) -l 



£ dt 



' T 27F 

ou i = — . 

UJ 

On definit egalement la moyenne spatiale de la densite d’energie : 



(£) = 



rx + A 



£ dx 



Dans le cas d’une onde progressive, ces deux valeurs sont egales et valent : 



<£> = 



vl 



2p 0 V* 



Intensite 

On appelle intensite de l’onde acoustique la puissance qui traverse, par unite de temps une 
surface unite perpendiculaire a la direction de propagation. 

Pour calculer l’intensite de l’onde calculons l’energie qui traverse pendant un intervalle de 
temps une surface S perpendiculaire a la direction de propagation. 



t t+dt 




:◄ ► 

V dt 

Flux de puissance 



Cette energie dE est egale a l’energie contenue dans un volume S V dt et elle egale a 

dE = £ SVdt 
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Niveau sonore 

On definit, le niveau sonore en decibels par 

N dB = 10 log 

Jo est une intensity de reference correspondant a Jo = 10 12 W ■ m~ 2 . 

Exemple : Pour une frequence / = 1kHz, le seuil d’audition est egal a 0 clB et le seuil de 
douleur est egal a 130 dB. Pour calculer l’intensite, l’amplitude de pression p 0 la vitesse de 
particules et le deplacement de particules, dans le cas on V = 330 m/s et Z c = 411 rayleighs, 
on pent utiliser les relations suivantes : 



I = I o l0 olNdB 
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8.5 Reflexion- Transmission 




Reflexion a un interface fluide-fluide 



Soit deux milieux Guides semi-infinis separes par un surface plane. Choisissons un repere 
orthonorme de telle sorte que le plan yOz coincide avec la surface de separation. Lorsque une 
onde acoustique provenant de — oo, se propageant. dans le premier dans la direction de l’axe des 
x arrive a la surface de separation, elle donne naissance a deux ondes 

- une onde reflechie qui se propage dans le premier milieu dans le sens des x decroissant. 

- une onde transmise qui se propage dans le second milieu dans le sens des x croissant. 
L’onde resultante dans le premier milieu (x < 0) est caracterisee par : 

Pi {x, t ) = pi (x, t ) + p R (x, t ) 

Pi(x,t) = Pi e j{u)t ~ kix) + P R e j ^ t+k ^ 

ui (x, t) = [Pi e j{ujt ~ kix) - P R e j{u)t+kix) ] 

Dans le deuxieme milieu, on a 

p 2 (X, t ) = pT (X, t ) 

p 2 (x,t) = P T e j{ - wt - k2x) 

u 2 {x,t) = ±-P T e?( ut - hax ') 

Z 2 

Les relations de continuite a l’interface s’ecrivent 

f pi(0,t) =p 2 (0,t) 

\ Ui (o ,t) = ii 2 (0 ,t) 

On en deduit 
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ou encore 




Pi + Pr — Pt 

(Pi - Pr) = P Pt 
A 2 



1 



P, 



R 



p P ' 



Pr 

)I\ 

Z\ Ft 

~Z~2~Pl 



On definit- 

- le coefficient de reflexion pour la pression 



Rp 



Pr 

Pi 



- le coefficient de transmission pour la pression 



T P = 

Les deux relations de continuity s’ecrivent alors 



Pt 

Pi 



f 1 + Rp — Tp 

l 1 -Rp = ¥Tp 

On en deduit les coefficients de reflexion et de transmission 



Z 2 — Z\ 

Z 2 + Z\ 

2 Z 2 

Z 2 + Z\ 

En tenant compte des relations p, = Z\ iii, pp = — Z\ v /,* et pp = Z 2 ut- on pent calculer les 
coefficients de reflexion etde transmission pour la vitesse de particules et pour le deplacement 
de particules : 



R P = 
T P = 



Rp 

Tu 



Ru — 



Z\ — Z 2 

Z\ + Z2 



T v = 



2Z\ 

Z\ + Z2 



En tenant compte des relations /, = Pf /2Zi 1 Ir = Pf ; j2Z\ et Ij = Pt/2Z 2 , on peut calculer 
les coefficients de reflexion et de transmission pour l’intensite acoustique : 
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or 

CXT 



ir = \ z!-z 2 y 

li _Z\ + Z 2 
It AZ\Z 2 
li [Z\ Z 2 
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Chapitre 9 



Propagation d’une onde electrique 
dans une ligne coaxiale 

9.1 Introduction 

Lors de l’etude des circuits electriques, en regime sinusoidal force, on considere habituelle- 
ment que dans une branche quelconque du circuit, le courant electrique est le meme en tout 
point et ne depend que du temps. Cette hypothese n’est valable que si la longueur d’onde 
associee au courant sinusoidal est grande devant les dimensions du circuit ; par exemple pour 
une frequence / = 50Hz, la longueur d’onde correspondante est A = c/ f ~ 6000/cm, ou c est la 
vitesse de propagation des ondes electromagnetiques dans le vide (c = 3 x 10 8 m ■ s _1 ). Cette 
hypothese revient de fagon equivalente a negliger le temps de propagation de l’onde electroma- 
gnetique (propagation de vitesse infinie). 

Si A est du meme ordre de grandeur des dimensions du circuit, cette approximation n’est 
plus valable ; par exemple si / ~ GHz, A ~ 30cm qui est alors comparable aux dimensions des 
circuits ou des fils couramment utilises dans les reseaux electriques a l’echelle du laboratoire. 
Dans ce cas, il faut tenir compte de la vitesse finie des ondes electromagnetiques. 



9.2 Equation de propagation 

Un cable possede une capacite par unite de longueur C\ et une autoinductance par unite 
de longueur L/ . Nous supposons que ce cable, que nous appelons ligne, est rectiligne et que la 
position de chaque point est reperee par son abscisse x. Soit un element de cable situe en x et 
de faible longueur dx. II possede une capacite Ci dx et une inductance L/ dx. Cet element de 
cable est represente, a un instant t, par le schema ci-dessous. 
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V'O ► 


—m — 


— ► — o- 4 


i(x) 


L /dx 


i(x+dx) 


v(x) 


C/dx — 


— 


■ -o 4 




r 

o ■■■■ 



v(x+dx) 



dx 



x 



H ► 

x+dx 



cellule elementaire 



1. L’ecriture des lois de Kirchhoff permet d’etablir les relations suivantes : 

dv ( x , t ) di ( x , t) 

dx 1 dt 

di {x, t) _ dv (a;, f) 
da; 1 dt 

A partir de ces deux relations, on peut etablir les equations de propagation de la tension 
v et du courant i : 



Dans ces deux equations, V 
dans le cable . 



d 2 v 1 d 2 v 

dx 2 V 2 dt 2 

d 2 i 1 d 2 i 

dx 2 V 2 dt 2 

es t v it(;sse de propagation de l’onde electrique 



9.3 Solution de l’equation de propagation 



On peut verifier que les fonctions / (t — y) et g (t — sont deux solutions de l’equation 
de propagation. La solution de l’equation de propagation s’ecrira done generalement sous la 
forme : 




La fonction f (t — y) correspond a une onde progressive se propageant dans le sens des x 
croissant tandis que g (t + correspond a une onde se propageant dans le sens des x decrois- 
sant. 
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9.4 Onde Progressive sinusoidale 

9.4.1 Definition 



Lorsque g (t + y) =0, la tension electrique en un point d’abscisse x est donnee par : 

v{x,t) = f (t- 

Dans ce cas, on dit qu’un onde progressive se propage dans le cable. 

Lorsque la fonction f (t — y) est une fonction sinuoi'dale, l’onde progressive est dite sinu- 
soidale ou harmonique : 

v (x, t ) = Vo cos 

On peut encore ecrire 

v (x, t) = Vo cos [cut — kx] 

on k = y (Relation de dispersion). 

En utilisant la notation complexe : 

v (x,t) = Re [n 0 e j[ut ~ kx] ] = Re [v e jut ] 

ou 1’ amplitude complexe est definie par : 




x : 



De meme , on peut obtenir l’expression de l’intensite en fonction du temps t et de la position 

i (x, t) = io cos (c at — kx) 



i (x, t ) = Re [i e jut ] 

ou 1’ amplitude complexe du courant est definie par la relation : 




9.4.2 Impedance en un point 

L’impedance en un point est definie par le rapport : 

v 
i 

Dans le cas d’une onde progressive, l’impedance est un nombre reel egal a 

Z c = 

et ne depend pas de x. Cette impedance est appelee impedance caracteristique ou impedance 
iterative du cable. 




Z(x) 
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Chapitre 10 

Elements d’analyse vectorielle 



10.1 Champ scalaire - Champ vectoriel 

Soit un triedre orthonorme (e x ,e y ,e z ) et M un point de l’espace, de coordonnees (x. y. z) : 

> _ 

OM = x e x + y e y + z e z 

La fonction / (M) est dite fonction scalaire de point ou champ scalaire si : 

f(M ) = f(x,y,z) 

Le vecteur v(M) est dit- fonction vectorielle de point ou champ vectoriel si : 

v (M) = V x (x, y, z) e x + V y (. x , y, z) e y + V z (x, y, z) e z 



10.2 Gradient d’un champ scalaire 

>• 

Le gradient ( note grad ) est defini a partir d’une fonction scalaire de point et a pour 
composantes suivant e x , e y , et e z les derivees partielles de / (M) par rapport ax, y et z respec- 
tivement : 

7> . ... Of _ Of Of 

grad (f) = - e I + - e s + - e> 



10.3 Divergence d’un champ vectoriel 



La divergence (notee div ) n’est definie qu’a partir d’une fonction vectorielle v (. M ) de point 
et donne une fonction scalaire de point definie, en coordonnees cartesiennes par : 



div (v) 



Ovx 

Ox 



Ovy 

0y 



OVz 

Oz 
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10.4 Rotationnel d’un champ vectoriel 



Le rotationnel note ( rot ) d’un champ vectoriel donne une fonction vectorielle de point 
definie en coordonnees cartesiennes par : 



rot (■ v ) = 



dv z dv y 
dy dz 



+ 



dv r dv~ 



dz dx 



dv y _ C)Vx_ 
dx dy 



10.5 Laplacien scalaire 

Le laplacien scalaire d’une fonction scalaire de point ( note lap ou A) est par definition un 
champ scalaire defini par : 

r — 

lap (/) = A / = div grad (/) 

Dans un systeme de coordonnees cartesiennes, il s’ecrit : 

<«p(/) = a/ = A + A + A 

dx 2 dy 2 dz 2 



10.6 Laplacien vectoriel 



Le laplacien vectoriel (note lap ou A) d’un champ vectoriel v est un champ vectoriel defini 
par : 

lap (v) = A (v) = grad [div (tT)] — rot rot (v) 

Dans le cas d’un systeme de coordonnees cartesiennes, le laplacien vectoriel a pour compo- 
santes : 

d 2 v T d 2 v r d 2 v r 



= dx 2 



+ 



lap (v) < 



AVy dx 2 



dy 2 dz 2 
y , ^ v y 



d 2 v v d 2 v n , d 2 u 



dz 2 



Av ’ = dx 2 



dy 2 

d 2 v z . d f v z . d 2 v. 



+ 



dy 2 



+ 



dz 2 



10.7 Operateur nabla 

Pour ecrire de maniere plus compacte les operateurs vector iels precedemment definis, on 
introduit un vecteur symbolique appele operateur nabla et defini par : 



d 



d 



d 



^ ^ x q "b 

dx dy dz 



Les operateurs vect-oriels s’ecrivent parfois a l’aide de l’operateur nabla sous les formes 
respectives suivantes : 
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le gradient d’un champ scalaire / est note 



Of „ , df 



dx ' dy ~ y dz 



grad (/) = V/ = — e x + — e y + 
la divergence d’un champ vectoriel est notee 



df 



div (v) = V • v = 

le rotationnel d’un champ vectoriel est note 

dv z dv v 



dv r dv v dv 7 

" + ~ + 



dx dy dz 



rot (v) = V x v = 



e x + 



dy dz 

le laplacien scalaire d’un champ scalaire est note 



dv x dv z 
dz dx 



dv y _ dv^ 
dx dy 



lap (/) = A f = div \grad (/) = V • V/ = V (/) 



V 2 se lit ’’del de”. 



le laplacien vectoriel d’un champ vectoriel est note 
lap ( v ) = A (v) = grad [div (tT)] — rot \ rot (v) 



= VV • (-e) — V x 



V x v\ 



10.8 Theoreme de Stokes 



10.8.1 Circulation d’un champ vectoriel 



On definit la circulation d’un vecteur v le long d’un contour ( C) . par l’integrale curviligne : 

B 




La circulation de long d’un contour ferme est notee : 
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10.8.2 Flux d’un champ vectoriel 

On definit. le flux d’un vecteur v a travers une surface (S') par l’integrale double : 



0/(S) (v) = JJ v-ndS 

( 5 ) 




Lorsque la surface (S') est fermee, le vecteur unitaire n est dirige de l’interieur vers l’exterieur. 



10.8.3 Theoreme de Stockes 

La circulation d’un vecteur le long d’un contour ferme ( C ) limitant une surface (S') est egal 
au flux de son rotationnel a travers cett-e surface. 



C (v) = <t>/{S) (r<& (v)) 
v ■ dl = JJ rot (v) ■ n clS 

( 5 ) 

le vecteur unitaire n est oriente selon la convention du tire-bouchon de Maxwell. 



10.8.4 Theoreme de Gauss- Ostrogradski (ou theoreme de la diver- 
gence) 

Le flux d’un champ vectoriel a travers une surface fermee (S')est egal a l’integrale de sa 
divergence dans le volume (r)limite par la surface fermee (S) : 



v-ndS 



div ( v ) dr 



(S fermee) 



W 
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Chapitre 11 

Les equations de Maxwell dans le vide 



11.1 Le champ elect romagnetique 

Proprieties du champ elect rost at ique 

Dans les etats stationnaires, le champ electrique est appele champ electrost-atique. Le champ 
electrostatique E cree par une distribution de charges de densite p situee dans le vide, est a 
circulation conservative, c’est-a-dire qu’il satisfait les relations integrate et locale : 

j)E - ell = 0 
(D 

rot (e\ = 0 

ou (r) est un contour ferme quelconque oriente. 

II satisfait egalement satisfait les relations integrale et locale : 

/ / s ' dS ^ Illi dT 

(S fermee) (r) 

(e) = £- 

ou (S') est une surface fermee quelconque orientee vers l’exterieur et (r) est le volume 
interieur a (5 1 ). e 0 = 8.854187817 x 1CT 12 F m _1 est la permittivite du vide 

11.1.1 Champ electromoteur et vecteur densite de courant 

Le champ electromoteur 

Lorsqu’un courant electrique circule dans un conducteur, cela implique Fexistence d’une 
force motionnelle f m agissant sur les porteurs de charge q et Fon debnit le champ electromoteur 
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E rn par la relation 





Q 



La circulation de ce champ le long d’un contour ferine oriente (T) n’est pas conservative 
c’est-a-dire qu’elle est differente de zero. Par definition cet-te circulation est appelee la force 
electromotrice e relative au contour considere : 




(T) 



On peut bien entendu definir la f.e.m relative a un trongon AB oriente, non ferine : 




Le vecteur densite de courant 

Un courant electrique correspond a des charges electriques mobiles. On appelle vecteur 
densite de courant j , le vecteur tangent a la ligne de courant, et defini par 

j = pv 

ou p est la densite volumique de charges mobiles et v la vitesse d’entrainement de ces charges 
mobiles. Le module de ce vecteur represente la charge qui traverse par unite de temps, l’unite 
de surface perpendiculaire a la direction de deplacement des charges mobiles ; il s’exprime en 
A.m~ 2 . 

En regime stat-ionnaire, c’est-a-dire lorsque le vecteur densite de courant j est independant du 
temps, le flux de j est conservatif ce qui se traduit par les relations integrale et locale : 

j ■ dS = 0 

( S fermee) 

div (j'] = 0 




11.1.2 Le champ magnet ique 

Le champ magnetique B cree par une distribution de courants de densite j est a flux 
conservatif, c’est-a-dire qu’il satisfait les relations integrale et locale : 




div 




0 

0 
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ou (.S') est une surface fermee quelconque. 

Le champ magnetique B satisfait les relations integrate et locale : 




(r) ( 3 ) 

rot (u) = /i 0 J 



ou (r) est un contour ferine et oriente quelconque et (S) une surface quelconque s’appuyant 
sur (r) et orientee a partir de (r) par la regie dite du "t.ire-bouchon de Maxwell" ou "du 
bonhomme d’Ampere". /i 0 = 1.2566370614 x 10 G NA 2 est la permeabilite magnetique du 
vide. 

11.2 Le regime variable 

11.2.1 Le phenomene de propagation 

Considerons un ensemble constitue par des circuits parcourus par des courants et par des 
distributions de charge variant en fonction du temps ; cet ensemble pouvant et-re au repos ou 
en mouvement. Au voisinage de ces distributions regnent un champ electrique et un champ 
magnetique. Contrairement au cas stationnaire, ces champs ne sont pas synchrones avec les 
sources, c’est-a-dire qu’a un instant t donne, ces champs dependent des valeurs des sources a 
l’instant t — 9 ; nous exprimons ce fait en disant qu’il y a propagation a vitesse finie des champs 
a partir des sources qui leur donnent naissance et le retard 9 est d’autant plus grand que le 
point ou Ton desire connaitre les champs est eloigne des sources. 

11.2.2 Le phenomene d’induction 

Un circuit filiforme au repos et parcouru par un courant invariable n’ent.raine l’apparition 
d’aucune f.e.m ou d’aucun courant dans un autre circuit filiforme au repos. II n’en est pas de 
meme si le courant varie ou si les circuits en presence se deplacent. l’un par rapport a l’autre : la 
f.e.m ou le courant qui apparaissent sont dus au phenomene d’induction. Ce phenomene entraine 
l’apparition d’un char up electrique supplementaire (appele champ induit) ; ce qui conduit a 
modifier la propriety fondamentale du champ electrique. 

11.2.3 Le phenomene de capacite 

Un circuit comprenant un condensateur alimente par une source de tension variable en 
fonction du temps, est parcouru par un courant variable bien que la continuity electrique soit 
interrompue par l’espace ent.re les armatures du condensateur. Dans ce cas l’intensite du courant 
n’est plus conservee tout au long du circuit puisqu’elle est nulle dans l’espace entre les armatures. 
II n’est done plus possible d’appliquer le theoreme d’Ampere. Pour conserver la validite de ce 
dernier, nous serons amenes a introduire un courant "fictif" appele courant de deplacement. 
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11.3 L’induction electromagnetique 

11.3.1 Loi de Faraday-Lenz 

On peut induire une f.e.m dans un circuit filiforme ( C ) ferme en faisant varier le flux 
magnetique a travel's le circuit : c’est. le phenomene d’induction electromagnetique. 

Pendant un temps dt. la variation du flux magnetique total a travel's une surface quelconque 
s’appuyant sur le circuit ( C ) est dd ; la f.e.m induite e s’exprime a l’aide de la loi de Faraday : 

d(j) 

dt 

Cette loi, etablie experimentalement pour des variations relativement lentes du flux magne- 
tique en fonction du temps, est valable pour tout regime variable et elle sert de base a l’etude 
de Telectromagnetisme classique. 

Puisque une f.e.m apparait dans le circuit ( C ) et y fait circuler un courant ceci implique 
l’existence d’un champ electromoteur agissant sur les porteurs de charge du circuit ((7). Ce 
champ est appele champ electrique induit. 

11.3.2 Equation de Maxwell-Faraday 

Considerons un circuit ( C ) au repos soumis a un champ variable. Un champ electrique 
va prendre naissance dans tout l’espace ou existe un champ magnetique variable. Le champ 
electrique induit joue un role de champ electromoteur et la f.e.m apparaissant dans tout le 
circuit ( C ) peut s’ecrire : 



e — lE,-dl^- d 4 = -d 
J dt dt 

(C) 



B-dS 



C s ) 



ou (,S') est une surface orientee s’appuyant sur le contour oriente ( C ). 
En permutant les operateurs d’integration et de derivation on a : 



f 



(J) Ei ■ dl = 

{C) 



d B-dS 



dt 



as ,s 

-m' ds 



(S) 



car le circuit ( C ) et.ant. immobile, la surface (.S') Test aussi et dS est independant du temps. En 
appliquant le theoreme de Stokes nous pouvons ecrire : 

Ei ■ dl = (4) -dS = -Jj^-dS 

(C) (5) (5) 

Cette egalite etant. satisfaite quel que soit ( S ) s’appuyant sur ( C ), il en resulte : 

dB 
~dt 



rot ( Ei) = 
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Remarquons que s’il existe en plus du champ electrique induit un champ electrostatique Eg. 
Le champ total E = Ei + Eg satisfait encore la relation ci-dessus car- 



rot = 0 



et on a la relation de Maxwell-Faraday 

rot (e) 



dB 

~di 



En definitive, nous devons retenir de l’etude du phenomene d’induct.ion electromagnetique 
le resultat fondamental suivant : en chaque point de l’espace ou existe un champ magnetique 
variable nous devons associer un champ electrique induit variable a circulation non conservative 
c’est-a-dire ne derivant pas d’un potentiel. L’ensemble de ces deux champs constitue le champ 
electromagnetique. 



11.4 Le theoreme d’Ampere 

11.4.1 Equation de continuity 

Si on considere une surface fermee (S) entourant un volume (r) . Si p est la charge volumique 
et q la charge totale du volume (r) a l’instant t ; on a 



Q = 




p dr 



(d 



Pendant l’intervalle de temps dt , la variation de la charge totale est dq et on a : 

/ \ 

dq d fff 

di = Tt J J p dr 

\ M / 



mais 



j-dS represente la charge totale sortant de la surface (S) par unite de temps, 



( S fermee) 

done par suite de la conservation de la charge : 



et 



d_ 

dt 



dq 

dt 



( 



J • dS 




(S fermee) 

\ 



p dr 

V ' W / 



j • dS 



( 5 ) 
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En changeant l’ordre des operations par rapport a l’espace et par rapport au temps, on a : 

'dp 




dt 



clr = — 



3 • dS 



(r) ( S fermee) 

or, d’apres le theoreme d’Ostrogradsky : 



j • dS = 



( S fermee) 




div (j\ dr 



done 




(t) 



% + dtv ( j 



dr = 0 



Cette relation doit etre verifiee quel que soit le volume (r), il faut done que Ton ait : 

iw ( 7 ) + 1 = 0 

Cette equation dite de continuite traduit la conservation de la charge electrique et montre 
que le flux du vecteur densite de courant n’est plus conservat-if comme dans le cas des et-ats 
stat-ionnaires. 

Ce qui se traduit localement par 

div (e)= — 

V / £ 0 

p represente la densite volumique de charges . L’equat-ion de continuite peut alors s’ecrire : 



div ( j ) + ^ = div j + ^ {div {e 0 E) j = div | j + ^ I = 0 



On voit que la relation fondamentale de conservation du flux de la densite de courant sera 
conservee si on fl applique a une densite de courant total Jt egale a la somme de la densite de 
courant vrai j (appele courant de conduction) et d’une densite de courant fictif appele courant 
de deplacement et defini par : 



3d — 



d z 0 E ) 



dt 



d’ou : 



jT=j+jD=j 



et div (j T ) =0 
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11.4.2 Le theoreme d’Ampere 



Le theoreme d’Ampere peut etre generalise a condition de l’appliquer au courant total.La 
relation de Maxwell- Ampere qui en est la traduction s’ecrit : 

— > / x OE 

rot [Bj= /i 0 j +li 0 £ 0 — 

La relation integrate du theoreme d’Ampere generalise est : 



/ 




ho 





d 




11.5 En resume 



En dehors des discontinuity, les equations de Maxwell s’ecrivent : 
- Theoreme de Gauss pour E 
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Chapitre 12 

Propagation des ondes 

elect romagnetiques dans le vide 



12.1 Equations de propagation pour E et B 

En absence de charges et de courant, les equations de Maxwell s’ecrivent : 



div ( E ) = 0 



( 12 . 1 ) 



div [ B ) = 0 



rot ( E ) = — - 



— * ( 

rot = g 0 e 0 — 

Pour etablir l’equation relative au champ electrique E. il eliminer le champ magnetique B 
. Pour cela, calculons le rotationnel de chacun des membres de le loi de Fraday : 



rot f rot (£'))= rot — 



(12.2) 



en permutant l’ordre des derivations, on obtient : 



Sachant que 



— * ( — * E\ ^ ( — >• E\ d f 9E 

rot (rot B) = [rot B) = -- W „— 



rot rot ( E ) = gracl div ( E) — A E 



(12.3) 



(12.4) 



ou A est le laplacien vectoriel, on obtient P equation aux derivees partielles suivante : 
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d 

A E - grad ( div [ E ) ) = — { g 0 e 0 



dt 



dE 

~0t 



Comme en dehors des charges 



on obtient finalement : 



div [ E ) = 0 



^ d 2 E > 

AE — g 0 e 0 -^- = 0 



(12.5) 

( 12 . 6 ) 

(12.7) 



Pour etablir Pequation aux derivees partielles pour le champ magnetique B. calculons le 
rotationnel de chacun des membres du theoreme d’Ampere-Maxwell : 



rot (rot j j = grad (div — AB = rot I /j, 0 £q 



dE 

~dt 



mais 

div (^B^j = 0 

et en inversant- l’ordre des derivations : 

— > d r— > 

-A B = +g 0 £ 0 — rot f E 



or 



Done 



ou encore 



rot ( E i = 



dB 



-AB = -fi 0 e o 



dt 

d 2 B 

~dB 



s, q2 ^ 

AB - h q eq = 0 

On obtient la meme expression que pour le champ electrique. 
L’equation 



E 

B 



i = o 

c 2 dB v B 



( 12 . 8 ) 

(12.9) 

( 12 . 10 ) 

( 12 . 11 ) 

( 12 . 12 ) 

(12.13) 



(12.14) 



c = 



constitue Pequation de propagation du champ electromagnetique dans le vide, ou on a pose 
1 

V% £ ° 
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12.2 L’onde plane progressive sinusoidale 

12.2.1 Relation de dispersion 

L’onde plane progressive sinusoidale est definie, en notation complexe, par : 

E (r,t) = E 0 e^-^) (12.15) 

ou k est le vecteur d’onde donnant la direction de propagation de l’onde plane. 

En utilisant la defintion du laplacien vectoriel dans un systeme de coordonnees cartesiennes, 
on peut montrer que : 

A E = -k 2 E (12.16) 

L ’equation de propagation s’ecrit alors sous laforme 

-k 2 +^r E = 0 (12.17) 

c z 

L’onde plane progressive sinusoidale constitue une solution particuliere de l’equation d’onde 
seulement si la relation suivante, dite relation de dispersion, est satisfaite : 

k = ^ (12.18) 

12.2.2 Structure de l’onde uniforme plane 

L’onde plane progressive sinusoidale doit egalement satisfaire le theoreme de Gauss. On 
montre aisement que pour une onde plane progressive sinusoidale : 

cliv = 0 est equivalents a — i k ■ E = 0 (12.19) 

Soit- encore k ■ E = 0 ; ce qui revient a dire que le champ electrique E est perpendiculaire a la 
direction de propagation donnee par le vecteur d’onde k. Le champ electrique est dit transversal. 




y 
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12.3 Onde plane uniforme progressive et sinusoidale 



12.3.1 Onde de polarisation rectiligne 

Jusqu’ici nous avons considere un champ electrique E dont la variation en fonction de z 
et de t etait quelconque. Nous allons maintenant examiner le cas ou le champ E garde une 
direction constante (polarisation rectiligne) et varie sinusoi'dalement ; il s’ecrit alors : 



E 



Eq cos 



cu 




E 0 cos [cut — k ■ r 



ou 

E 0 = vecteur constant 

cu = pulsation de la fonction sinusoidale 

k = lu/c : constante appelee nombre d’onde 

k = k n : vecteur d’onde perpendiculaire au plan d’onde 

cut — k ■ f : phase instantanee ou plus simplement phase de la grandeur variable, 
c = - j== vitesse de propagation dans le vide . 

Les vecteur s ferment toujours un triedre direct et 

-> it x E k x E 

B = = 

c UJ 

Pour preciser cette onde, supposons qu’elle se propage suivant z'z d’ou 

E = E 0 cos (cut — kz ) 



On constate une double periodicite : 



Une periodicite temporelle : pour z donne le champ varie sinusoi'dalement en fonction du 
temps avec une periode 

T — — 



cu 



ou une frequence 




cu 

2n 



(/ est en hertz). 

-Une periodicite spat-iale : a un instant t donne le champ varie sinusoi'dalement en fonction 
de z avec une periode 

2tt _ A 
~~k~k 



(A est appelee la longueur d’onde dans le vide). On peut remarquer que la longueur d’onde 
A est egale a la distance parcourue par l’onde pendant une periode. 



http:/ / dj elouah.ifrance.com 






12.3 Onde plane uniforme progressive et sinusoidale 



94 



12.3.2 Onde de polarisation quelconque 

Dans le paragraphe precedent, nous avons suppose que le champ E (done B egalement) gar- 
dait- une direction const-ante. Dans le cas general, il n’en est pas t-oujours ainsi et les composantes 
du champ peuvent se mettre sous la forme : 



E x = E 0x cos (cut — kz — 4i) 

E y = E 0y cos (cut — kz — 44) 

0 i, 0 2 pouvant etre differentes. 

Etudions le comportement du champ E dans le plan z = 0. Les result-ats obt-enus se ret-rouvent 
avec un decalage t-emporel dans tout plan z = cte. Les composantes du champ s’ecrivent- : 



E x = E 0x cos (cut - 44) 

E y = E 0y cos (cut — 44) 

E z = 0 

et si Ton prend pour origine des temps un instant ou E x passe par sa valeur maximale on 
a : 



E x 

Eq x 



E 0 y 



cos (cut) 
cos (cut — 4) 



avec 4 — 02 0i • 

On peut, deja dire que l’extremite du vecteur decrit une courbe inscrite dans un rectangle 

E 

de cotes 2 E 0x et 2E 0y . D’autre part en developpant- l’expression de et en eliminant- le temps 
il vient : 



Eo. 



E, 



o y 



= cos (cut) cos (0) + sin (cut) sin (0) 



E n , 



E, 



0 y 



E X COS (0) + \ l 1 - ^ 



Ei 



Ox 



E, 



sm i 



v , 



r ee 
7 ^ x cos(0) 


2 




_E 0 y E 0x 




V-^Ox / 



snr 



E x 



Ei 



Ox 



En. 



Ei 



0 y 



E x E v 

2 f^f^ cos{ 



= sin 2 (0) 



Pour 0 quelconque, cet-t-e equation est celle d’une ellipse : on dit- que l’onde a une polarisation 
ellipt-ique ; pour 0 = mn(m = 0, 1, 2, ... ) l’ellipse degenere en une droit-e et l’onde est dit-e a 
rectiligne. Enfin si E 0x = E 0y et si 0 = (2 1 + l)n/2 l’onde est dit-e a polarisation circulaire. 
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12.3.3 Notation complexe 

Considerons une onde plane sinusoidale polarisee rectilignement se propageant suivant une 
direction n et ayant pour vecteur d’onde k = k n . 




Le support du vecteur champ electrique E de direction constants a pour vecteur unitaire u 
et nous pouvons ecrire : 

E = E 0 cos (cut — k ■ r + ip'j 

sous la forme 

E = E 0 cos (cut — k • r + i/A u 



Avec Eq = cte =amplitude du champ E. 

On sait qu’une grandeur sinusoidale peut etre representee par un nombre complexe et qu’a 
chaque instant cette grandeur est obtenue en prenant la partie reelle du complexe qui la repre- 
sente. Ainsi dans le cas present, nous pouvons representer la mesure 

E = E 0 cos (cut — k ■ r + ib) 



du vecteur champ electrique sur son support par le nombre complexe 

E* = Eq exp (cut — k ■ f + 
et nous avons bien a chaque instant 

E = Re (E*) 

Ce complexe appele la composante complexe de E sur son support peut aussi bien s’ecrire : 

E* = En = Jj e jut 



E = Eq 
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est Pamplitude complexe de la composante de E sur a . Finalement nous caracterisons une 
onde plane polarisee rectilignement par sa representation complexe : 



E* = E 0 e j ^ 
B* = B 0 e ji> 



12.4 Energie electromagnetique : vecteur de Poynting 

La propagation de l’energie se manifeste experimentalement dans de nombreux cas : 

'f* On peut ressentir son effet si Pon s’expose aux rayons solaires ou au rayonnement d’une 
source chaude ; 

- De meme tout emetteur radio expedie de Penergie a travers Pespace, une infime partie de 
cet-te derniere etant- captee par votre recepteur radio. 

Nous allons essayer de relier localement cette energie qui se propage, au champ electro- 
magnetique qui la transports. Nous supposerons le milieu de propagation parfait, c’est a dire 
homogene, isotrope et lineaire. 



12.4.1 Onde de forme spatiale et temporelle quelconques 

Nous admettrons que les densites d’energie electrique et magnetique calculees en regime 
stationnaire sont toujours valables en regime variable ; la densite d’energie electromagnetique 
w en un point quelconque du milieu parcouru par une onde electromagnetique est done a chaque 
instant : 

w = densite d’energie electrique + densite d’energie magnetique 





Considerons dans le milieu, un volume r limits par une surface (S). L’energie electroma- 
gnetique qu’il contient est a chaque instant : 



W = 




Pendant un temps dt l’accroissement d’energie dans (r) sera dW et la puissance inst-antanee 
p' acquise par ce volume sera 

, dW 



P = 



dt 




dw 

~dt dT 



(■ T ) 



On a : 



\ dB 


^ 


B 


' “ ~~dt 


et rot 


ho 



— £o' 



dE 

~dt 
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done 



dw 

~dt 



= E ■ 



rot 



B 

do 



B 

do 



■ rot E 



D’apres une relation de transformation, on a : 

div 



done 



et 



Ex — 
do 



= — • rot (e) - E ■ rot ( — 



do 



do 



diu / - B 

—— = —div \ E x — 

dt \ /i 0 



P = 




(t) 



B 

div | Ex — | c/r 
do 



La puissance electromagnetique instantanee perdue par le volume (r) est : 



~P = 




(t) 



B 

div | fix — I c/r 
do 



Elle represente la puissance electromagnetique qui sort du volume (r), e’est a dire la puis- 
sance moyenne p rayonnee par ce volume. 



P = 




(t) 



B 

div | Ex — | c/r 
do 



D’apres la formule d’Ostrogradsky, on peut ecrire : 



Le vecteur 




• dS = 




P-dS 



( 5 ) 



R 



E 



B 
x — 
do 



est appele le vecteur de Poynting. Sa direction donne en chaque point, la direction d’ecoule- 
ment de l’energie et son flux a travel's une surface est egal a la puissance electromagnetique 
instantanee rayonnee par cette surface. Les courbes tangentes en chaque point au vecteur de 
Poynting peuvent etre considerees comme des trajectoires de l’energie ; on les appelle les ‘’rayons 
electromagnetiques” . 
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12.4.2 Onde plane progressive et uniforme sinusoi'dale 

Puisque ^E, B, fi\ forment un triedre trirectangle direct le vecteur R a meme direction et 

sens que k c’est a dire que l’energie s’ecoule dans le sens de propagation (ce result-at n’est pas 
general ; en effet dans un milieu anisotrope par exemple R et k ne sont pas colineaires) . 

La puissance instantanee p„ traversant une surface unitaire (S) perpendiculaire a la direction 
de propagation est 



Pu= j j R ■ ds = 

( 5 ) ( 5 ) 



R 



dS = 



R 



clS = 



R 



S 



( 5 ) 



La puissance moyenne traversant (S) est alors 



or B _L E, et 



B 



T 



— \/do £ o 



1 f T 




1 


f T 


- B 


dt = — / 


R 


dt = — 




Ex- 


T Jo 




T 


Jo 


do 



dt 



E 



, d’ou 



(Pu} = 



1 

T 



/ E\^dt 
I o V / J o 



Si l’onde est polarisee rectilignement alors 

E = E 0 cos [cut — k ■ r 



(p«) 



TP f \ ^ cos2 — k-r\ dt 

T Jo V do V J 



(Pu) 




f T 1 


r 


2 (ojt — k ■ J ) 


/ h 


1 + COS 


Jo 2 


- 


. V / J . 



dt 






%1L 

z n 



ou E e f f = valeur efficace de E. 

Le flux d’energie traversant par unite de temps l’unite de surface perpendiculaire a la di- 
rection de propagation est une constants dependant du milieu et proportionnelle au carre de la 
valeur efficace du champ electrique. 
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Chapitre 13 



Reflexion et transmission des ondes 
elect r omagnet iques 



Nous allons dans ce chapitre utiliser les equations de maxwell pour etudier la propaga- 
tion d’une onde d’une onde electromagnetique dans un milieu isolant puis sa reflexion et sa 
tranmission a un interface. 



13.1 Equations de Maxwell dans les milieux par fait s 

Nous supposerons que l’onde se propage dans un milieu illimite satisfaisant les conditions 
suivantes ( milieu parfait) : 

•Le milieu est considere comme continu ; ce qui est legitime pour les ondes habituellement ren- 
contrees (hertziennes et lumineuses) dont- les longueurs d’onde sont, en general, tres superieures 
aux distances inter-atomiques. 

•Le milieu est homogene et isotrope 

•Le milieu peut et-re caracterise par une permittivite £ et un permabilite /i et une conductivity 
7 constantes egales a celles definies en regime stationnaire (milieu lineaire). 

Ces hypothese simplificatrices permettent de traiter un grand nombre de problemes ; il faut 
toutefois signaler qu’elles ne tiennent pas compte des phenomenes d’absorption, d’hysteresis et 
de dispersion ( £ , // et 7 peuvent alors etre complexes et dependre de la frequence). 

Dans le cas des milieux homogenes, lineaires et isotopes, les equations de Maxwell de- 
viennent : 

- Theoreme de Gauss pour E 



Forme locale 


Forme integrale 


div (e) = Plibres 


J 


f J E-dS = JJJ Plibres dr 

( S fermee) (r) 



- Theoreme de Gauss pour B 



http:/ / dj elouah.ifrance.com 



99 






13.2 Propagation dans les milieux dielectriques 



100 



Forme locale 


Forme integrale 


cliv (b) 


1 =0 


J 


f J B ■ dS = 0 

( S fermee) 



- Loi de Faraday 



Forme locale 


Forme integrale 


rat (^E^j 


i)B 
1 “ ~~dt 


J 

(C 


f E-dl = -j t JJz-dS 

) (S) 



- Theoreme d’Ampere-Maxwell 



Forme locale 


Forme integrale 


rot (^B^j 


dE 

1 = h 3 + VZ-gjT 


j 

(r 


Ib-ciT= JJ | 

) (5) 


( , BE) 

^3+l^j 


I • (IS 



ou 

- £ = e 0 £ r est la permittivite absolue du milieu dielectrique 

- £ r est la permittivite relative du milieu dielectrique. 

- £o = 8.854187817 x 10” 12 F m” 1 est la permittivite du vide 
/i = Podr es t la permeabilite absolue du milieu 

- /i r est la permeabilite relative du milieu. 

- /i 0 = 1.2566370614 x 10” 6 NA” 2 est la permeabilite magnetique du vide 



En tout point d’un conducteur, il existe une relation entre le vecteur densite de courant j et 
le champ electrique total E ( somme du champ electrostatique Eg et du champ electromoteur 
E rn ) dite relation d’Ohm Kirchhoff : 



j = 7 (Es + Em') 

3 = 7 E 

ou 7 est la conductivity du milieu conducteur. 

Cas particuliers : 

- Isolants : 7 = 0 

- Conducteurs parfaits : 7 — > 00 

13.2 Propagation dans les milieux dielectriques 

Les milieux dielectriques sont des mieux isolants. Leur conductivity est extremement faible, 
de l’ordre de 10 20 a 10 12 celle d’un conducteur metallique etant de l’ordre de 

10 7 (4) • m” 1 ), a temperature ambiante. II est done tout a fait raisonnable de prendre pour 7 
la valeur 7 = 0. Par ailleurs dans un tel milieu, Pi ibre = 0. 

Les equations de Maxwell se simplifent alors en : 
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div 


\ s ) 


= 0 


div 


(5) 


= 0 





( e) 


dB 


rot 




V ) 


dt 


> 


(5) 


dE 


rot 


= ^ 



En utilisant la meme demarche que dans le chapitre precedent, on peut montrer que le 
champ electrique et le champ magnetique satisfont les equations de propagation suivantes : 



A E 



A B 



1 d 2 E 
v^ljt 2 
1 d 2 B 



V 2 dt 2 

ou la vitesse de propagation de l’onde est : 



0 

0 



yj . yj f^O^Oy/ fJ> r E r Tl> 

n = v /// r c r est l’indice de refraction (ou indice optique) du milieu. Dans les milieux reels n 
est constant pour les grandes longueurs d’onde, t.andis que pour les hautes frequences il faut, 
faire intervenir le phenomene de dispersion qui entraine une dependance de n avec la frequence. 
Dans la plupart des dielectriques //,, = 1. d’ou n = \JE r . 

On peut egalment. montrer que l’impedance caracteristique d’un tel milieu peut s’ecrire : 



Z = 

ou Z 0 est l’mpedance caracteristique du vide. 



Zo 

n 



13.3 Relations de passage 

Composante tangentielle et composante normale de E 

A la traversee d’une surface (E) separant deux milieux et portant des charges vraies avec 
une densite supercielle a, les relations locales s’ecrivent : 



Et 2 — Eti 
£2 En 2 £1 Em 



0 

a 



ou E t est la composante de E dans le plan tangent a (E) en M , tandis que E N mesure de 
la composante de E suivant la normale n en M et orientee du milieu (1) vers le milieu (2). 
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Composante tangentielle et composante normale de B 

A la traversee d’une surface (S) separant deux milieux et parcourue par des courants vrais 
de densite superficielle js , les relations locales B deviennent : 

Bn2 = B n \ 

Ih Bt2 ~ Pi Bti = js x n 

ou B n est la composante de B suivant la normale au point M considere orientee du milieu 
( 1 ) vers le milieu ( 2 );tandis que B T est la composante de B dans le plan tangent a (£) en M. 

13.4 Lois de Snell-Descartes 

Considerons deux milieux 1 et 2 , parfaits, isolants, non magnetiques et semi-infinis, carac- 
terises respectivement par £\^i l = /j 0 , = 0 et par £ 2 ,Ai 2 = l L o> 72 = 0 ? et separes par une 

interface plane. On constate experimentalement qu’une onde plane incidente arrivant a cet-te 
interface donne generalement naissance a deux ondes, une onde plane reflechie et une onde 
plane transmise ou refractee ( la reflexion entre deux milieux parfaits et isolants est souvent 
appelee reflexion vit reuse ) 

Nous nous proposons de determiner entierement les ondes reflechie et transmise en appli- 
quant- les equations de propagation et en tenant compte des conditions aux limites. Quel que 
soit l’etat de la polarisation de l’onde incidente, il est toujours possible de la decomposer en 
deux ondes a polarisation rectiligne perpendiculaires entre elles en projetant le vecteur inci- 
dent sur deux axes perpendiculaires ; nous etudierons done d’abord le cas d’une onde incidente 
polarisee rectilignement. 

Definissons un triedre de reference par une origine O appurtenant, au plan separant les deux 
milieux, l’axe Oz etant perpendiculaire a ce plan et l’axe Oy et-ant tel que le plan yOz contienne 
le vecteur k\ de l’onde incidente. En appelant u 1 , u \ el - u 2 les vecteurs unitaires des supports 
des vecteurs champs electriques des differentes ondes, nous pouvons representer ces dernieres 
en notation complexe par 

- Onde incidente : Ei = E 0 i Tt 1 k i- r 1) 

- Onde reflechie : E[ = E’ m e u(^d-^'r^>) 

- Onde transmise : E 2 = E 02 it 2 e J ( a;2t_i:2 ' r2 ) 

Dans ces expressions E m , E ' 01 et E 02 sont les amplitudes complexes des differents champs 
electriques. En un point quelconque M 0 du plan xOy tel que OM 0 = r o,les ondes sont 
representees par 

- Onde incidente : Ei = E m iti k 1- r °) 

- Onde reflechie : E[ = E’ m lt[ 

- Onde transmise : E 2 = E 02 it 2 e J ( a;2t_i:2 ' ro ) 
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Reflexion-transmission 



La continuity de la composante tangentielle du champ entraine une relation entre les com- 
posantes tangentielles de la forme 

~a 1 E 01 A Ult - tl - r °) + ~a[ EoxA^-^ 0 ) = ~a 2 

Pour que cette relation soit satisfaite a chaque instant en tout point M 0 de la surface de 
separation, il faut que 



uj\t — k i • r o = co\t — k [ ■ r o = co 2 1 — k 2 ■ r 



ce qui implique 



uj\ = LOi = cu 2 et que k i • r 0 = k \ ■ r 0 = k 2 ■ r 0 
La derniere relation conduit a 

( k[ - k i) • rqj = 0 
(k 2 — k -~r o = 0 

qui seront satisfaites quel que soit le point M 0 si ((k [ — k ^ et (^k 2 — k sont per- 

pendiculaires au plan xOy, c’est a dire si k ±, k \ et k 2 sont coplanaires : ces trois vecteurs 
appartiennent au plan d’incidence ( k i,’e > 3 V 
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En prenant r 0 = e 2 , on obtient 

k i • ~e 2 => K sin (i\ ) = k\ sin (q) 
k i ■ ~e 2 => k2 sin (i 2 ) = k\ sin (q) 

u> 0J cu 

or ki = — , k[ = — et k 2 = ; d ou 

Vl Vl V 2 

sin (i^) = sin (q) => i\ = i\ 

sin (i 2 ) sin(?i) 

V^2 = ^ 

En multipliant les deux membres de la derniere relation par c (vitesse de propagation dans 
le vide) et en faisant- intervenir les indices des milieux, on obtient 

n\ sin (i 2 ) = n 2 sin (ii) 

En definitive nous avons montre que : 

- Les ondes reflechie et transmise ont meme frequence que l’onde incidente 

- Les vecteurs d’onde reflechi et transmis sont dans le plan d’incidence (plan et la normale 
a la surface de separation au point d’incidence) 

- L’angle de reflexion z^est egal a bangle d’incidence i\ 

- Les angles de transmission i 2 et d’incidence i{ sont lies par rq sin (z 2 ) = n 2 sin (q) (loi de 
la refraction). 

Ces resultats representent les lois de Snell-Descartes pour la reflexion et la refraction. 

13.5 Formules de Fresnel 

Nous allons etudier successivement le cas d’une onde dont. le champ electrique est polarise 
rectilignement dans le plan d’incidence puis le cas d’une onde dont le champ electrique est 
polarise perpendiculairement au plan d’incidence 

13.5.1 Champ electrique dans le plan d’incidence 

Le triedre de reference est choisi comme precedemment et les champs positifs des champs 
electrique et magnetique sont indiques sur la figure (Pour satisfaire les raisons de continuite 
et de symetrie, il faut supposer que les ondes reflechie et transmise ont leur vecteur champ 
electrique polarise dans le plan d’incidence). 

Sachant que les differents vecteurs unit.aires u \. u[ e t u 2 sont dans le plan d’incidence, les 
differentes ondes ont pour representation complexe : 

- Onde incidente : Ei = E 0 i ~u 1 e K UJlt ~^ 1 " r x ) 

- Onde reflechie : E[ = E ’ 01 ^'("i^i^i) 



k i • "e’ 2 = 
k 2 ■ ~e 2 = 
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- Onde transmise : E 2 = E 02 it 2 e j ( U2t k2 ' r 2 ) 




Champ electrique parallele au plan d’incidence 



Dans le milieu l l’onde resultante est la somme de l’onde incidente et de l’onde reflechie. Les 
continuities des composantes tangentielles du champ electrique et de l’excitation magnetique a 
la traversee de la surface separant les deux milieux permet, si Ton tient compt-e de // , = /i 2 = /j 0 . 
de B 0 = Eq/V et de la condition k 1 ■ r 0 = k [ ■ r 0 = k 2 • d’obtenir les relations 

E 0 i cos (*i) + E' m cos(ii) = E 02 cos (i 2 ) 

Eqi _ E ' 01 _ Eq 2 

“ vT 

(le signe moins provient du sens de B\ oppose a celui de Ox ; ce sens a ete choisi ainsi pour 
respecter le fait que ^E, B , k \ est toujours direct). 

A partir des relations ci-dessus, on obtient E[ n et E (y2 en fonction de E m : 



or 



-E'01 



cos (A) , p, 

cos(A) + 01 



E, 



02 



Vy 

Vo 



E' 



01 



E'02 



cos (i 2 ) 
cos (ii) 



E'01 



V\ _ n 2 _ sin (A) 
V 2 n\ sin (i 2 ) 
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d’ou les for mules de Fresnel 

2 cos (*i) sin (z 2 ) 2 cos (zi) sin (z 2 ) 

sin (ii) cos (zi) + sin (z 2 ) cos (z 2 ) sin (zi + z 2 ) cos (zi — z 2 ) 

sin (zi) cos (zi) — sin (z 2 ) cos (z 2 ) tan (zi — z 2 ) 

sin (zx) cos (zi) + sin (z 2 ) cos (z 2 ) tan (zi + z 2 ) 

Nous notons [ ]y pour indiquer que le champ electrique est dans le plan d’incidence. 

13.5.2 Champ electrique perpendiculaire au plan d’incidence : 

Les dispositions et les sens positifs pour les champs sont indiques sur la figure cidessous 





Champ electrique perpendiculaire au plan d’incidence 

si Ton tient compte de /j , = /j 2 = /i 0 . de B 0 = E 0 /V et de la condition k \-~r q = k' 1 -J^o = 
k 2 • r o les conditions de continuite a la traversee de la surface separant- les deux milieux 
permettent d’obtenir les relations 

E 01 + E' 01 = E 02 

Eoi / . x . -E'01 / • \ E 02 * . . 

— — cos(zi) + — cos(zi) = — — cos(z 2 ) 
vi V\ V2 

En tenant compte 

Y± _ n2 
V 2 n x 

on obtient finalement. les formules de Fresnel 



sin (zi) 
sin (z 2 ) 
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T ± 



R± 



Eq2 



Eqi 




Eqi 



2 sin (i 2 ) cos (L) 
sin (*i + i 2 ) 
sin [ii - i 2 ) 
sin (ii + i 2 ) 



Nous notons [] ± pour indiquer que le champ est perpendiculaire au plan d’incidence 



13.5.3 Discussion des resultats 

Dans les formules de Fresnel, les coefficients de reflexion et de transmission sont des nombres 
reels ; il en resulte que le dephasage introduit par la reflexion ou la transmission aura toujours 
une valeur nulle ou egale a n. 



Incidence normale : 



Nous obtiendrons les resulats de l’incidence normale en faisant- tendre ?'/ vers zero et en 
tenant compte alors de rq / 1 = n 2 i 2 . Les relations de Fresnel s’ecrivent alors 



T\ 

Ri 



T\ 



R± 



2 Ti l 
rii + n 2 
ni - n 2 
rii + n 2 



- Les ondes reflechie et transmise sont independantes de la polarisation de l’onde incidente 

- L’onde transmise est toujours en phase avec l’onde incidente 

- La composante electrique est dephasee de n a la reflexion si rq < n 2 ; elle n’est pas 
dephasee si rii > n 2 . 

Interessons nous maintenant- au comport-ement energetique de l’onde a l’interface plane ; ce 
comportement sera caracterise par les deux coefficients suivants 

T . puissance reflechie 

- Le pouvoir reflecteur an = 

puissance incidente 

. . puissance transmise 

- Le facteur de transmission = 

puissance incidente 

Les relations donnant, les intensites, c’est a dire les puissances par unite de surface, trans- 
portees par les difierentes ondes sont donnees par 



_ ffpi»-l . rf 

2 c/i ’ 1 



E&n x 

2 °Hi 



et I 2 



Eq 2 ti 2 
2 cn 2 



avec Hn 



l J 2 — L L o 



d’ou 



OiR 



a x 



'Ki 

_Eqi 

Eq 2 

Eqi 



1 2 



1 2 



ni — n 2 1 2 
_rq + n 2 _ 
4nin 2 
(ni + n 2 ) 2 
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Remarque : La conservation de la puissance implique or + cer = 1 ce que Ton peut verifier 
a partir des deux dernieres relations. D’autre part, le pouvoir reflecteur et le facteur de trans- 
mission ne dependent pas du sens de propagation de l’onde incidents (milieu 1 — > milieu 2 ou 
milieu 2 — > milieu 1). Pour les ondes lumineuse, dans le cas ou les milieux sont Pair et le verre 
(rq = /, ri 2 = 1.5) on a or = 0.04 et = 0.96. 

Incidence oblique 

Nous allons discuter le comport-ement de l’onde en fonction de son angle d’incidence q 
(0<*i < |) dans le cas particulier ou rq < n 2 : 

- Onde transmise : Puisque sin(i 2 ) = (^sin(q) on a * 2 < i\ et il existera t-oujours une onde 
transmise. II en resulte que le champ electrique de l’onde transmise, aussi bien dans le 
cas de polarisation dans le plan d’incidence (note ||) que de polarisation perpendiculaire 
au plan d’incidence (note _L), reste en phase avec celui de l’onde incidente sur la surface 
de separation 

- Onde reflechie : Pour une onde polarisee perpendiculairement au plan d’incidence (_!_) le 
champ electrique reflechi est t-oujours en opposition de phase avec celui de l’onde incidente. 
Pour une onde polarisee dans le plan d’incidence (||) et sous incidence normale, le champ 
electrique reflechi est en opposition de phase avec le champ incident. Lorsque Tangle 
d’incidence at-teint la valeur %\b definie par %\b + *2 = § , le champ electrique reflechi 
s’annule. Cet angle i\n est appele angle de Brewster et en tenant compte de n/ sin(q) = 
n 2 sin(z 2 ), on t-rouve qu’il a pour valeur i\B = atan(n 2 /?q). Pour un angle d’incidence 
proche de | (incidence rasant-e) ; le champ electrique reflechi a meme amplitude que le 
champ incident et ils sont en opposition de phase. 

Si l’onde lumineuse arrive directement- d’une source (le soleil par exemple) sur une surface 
vitreuse, elle est polarisee elliptiquement- et on peut la considerer comme la superposition 
de deux ondes polarisees rectilignement Tune dans le plan d’incidence T autre perpendicu- 
lairement au plan d’incidence. Suivant- Tangle d’incidence, la lumiere reflechie sera plus ou 
moins polarisees rectilignement cette polarisation etant t-ot-ale pour Tangle de Brewster 
la vibration lumineuse est- alors rectiligne et parallele a la surface reflechissant-e. 

- Pouvoirs reflect-eurs : Les pouvoirs reflect-eurs (int-ensit-e reflechie sur int-ensit-e incidente) 
ont respectivement pour valeur 



«i?|| 

®R± 



tan 2 (ii — i 2 ) 
tan 2 (ii + i 2 ) 
sin 2 (ii - i 2 ) 
sin 2 (ii + i 2 ) 



La figure dessous donne pour une onde lumineuse et dans le cas air- verre, la variation des 
pouvoirs reflecteurs en fonct-ion de Tangle d’incidence (i\B = 56,3°) . 
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Angle d'incidence 



Variation du pouvoir reflecteur en fonction de l’angle d’incidence 



13.6 Reflexion sur un conducteur parfait 



Considerons une onde plane incidents uniforme sinusoi'dale polarisee rectilignement ( E / || 
Ox) se propageant. dans le vide ( ou Fair) suivant zO et arrivant sous incidence normale, a la 
surface plane d’un conducteur de conductivity infinie (conducteur parfait). Des considerations 

energetiques montrent qu’il n’y pas d’onde transmise i^E 2 = 0 etB 2 = 0 j et que les seuls 
courants vrais pouvant etre induits par cette onde le sont sur la surface du metal. Les raisons 
de symetrie impliquent que les directions de polarisation des vecteurs incident et reflechi sont 
identiques. 




Reflexion metallique en incidence normale 
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Les champs electrique et magnetique des ondes incidents et reflechie sont 

~E ! = E 01 ~e 1 -,l3 l = e j{ “ t+kiz) ~e 2 

c 

E[ = E' m e j{wt - kiz) ~e 1 -,^ 1 = ^ e j{ut - kiz) ~e 2 

car ~u i = ~u\ = et , ~v i = ~v\ =e2; k i = — ki ~e 3 ; ~r i = ~r\ = z ~e $ et B 0 = — . 

c 

La continuity de la composante tangentielle du champ electrique permet d’ecrire a la tra- 
versee de la surface de separation (z = O ) 

{Eoi + E{ n } z=Q = 0 d’ou Eqi = —E ' 0 1 



Le champ electrique reflechi a meme amplitude que le champ electrique incident et il est 
dephase de n par rapport a ce dernier. Pour satisfaire a la directivite du triedre (^E , B , k j le 

champ magnetique est reflechi suivant B 1 c’est a dire que sa reflexion s’effectue sans changement 
de phase. En definitive les ondes incidente et reflechie s’ecrivent 



Et = E m e j{ “ t+kiz) ~et ; ~Bt = e j(ujt+kiz) ~e 2 

c 

E\ = —E 0 i e j{u}t ~ kiz) ~e x ; ~B X = - — e j{ut ~ kiz) ~e 2 

c 

Dans le vide, l’onde resultante est la somme de l’onde incidente et de l’onde reflechie et ses 
vecteurs champs ont pour valeur 

E = Et + E; = E 01 ( e +jklZ - e ~ jklZ ) e jut ~et 
E = 2E 01 sm(ktz) e^+i) ~et 



~B = Bt + B't = (e +jkiz + e~ jkiz ) e juit ~e t 

c v ' 

~B = 2 ^ cos (ktz)e j ^ t+n) ~et 

En revenant. a la notation reelle 

E = 2E01 sin (ktz) cos (cut + — j ~e 1 
^ 

B = 2 ^— cos (ktz) cos (cut + tt) et 

Ces relations montrent que les champs ne se propagent plus mais qu’ils oscillent sinusoi'dale- 

ment en fonction du temps avec une amplitude qui est fonction de la distance z, leur dephasage 
7 r 

et.ant de — on dit que 1 ’ onde est stat.ionnaire. 

Zj 
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L’amplitude du champ electrique est nulle pour k\Z = ‘=fz = l ir(l = 0,1,2, ) , c’est a 

dire pour z = 

L’amplitude du champ magnetique est uulle pour k\z = ^ z = (l + |) n (1 = 0, 1, 2, ) , 

c’est a dire pour z = (l + |) 

Les points ou l’amplit-ude est nulle sont appeles les noeuds et les points ou l’amplitude est 
maximale sont appeles les ventres. 
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